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ABSTRACT 


The variability of the marine atmospheric boundary layer (MABL) and its effect on 
the electromagnetic (EM) refractive structure around the Greenland Sea marginal ice 
zone Were examined. Rawinsonde profiles and surface observations collected from three 
ships during MIZEX-87 (20 March-11 April) served as the data set for the examination. 
A program was developed to calculate the refractivity at each vertical level of the 
rawinsonde profiles. The program also identified the levels at Which trapping, super- 
refraction and subrefraction occurred. 

Temporal studies showed that a higher incidence of anomalous refractive layers oc- 
curred during periods when the region was under the influence of high pressure. More 
than 50% of the time, trapping and super-refractive lavers were attributed to the devel- 
opment of a capping inversion just above the MABL during these periods. 

Spatial studies showed that the refractive structure varied relative to distance from 
the ice edge as did the depth of the MABL. An upward slope in refractive layer heights 
Was observed from the ice toward the open water. Significant spatial inhomogeneity was 
observed over horizontal ranges of less than 100 km. This was attributed both to the 
large-scale synoptic forcing affecting the region and to variations in the surface fluxes 
of heat and moisture over the ice and over the water. A range-dependent ray trace 
model developed at the Naval Ocean Systems Center was used to show how the ray 
paths of EM waves vary with a changing refractive structure. 
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I. INTRODUCTION 


A. GENERAL 

The Arctic region depicted in Fig. J is entering a period of greatly increased eco- 
nomic, social, strategic and political change. The major force behind the changes taking 
place is resource development activity. Although the present scale of this activity is not 
inconsequential, it is small in comparison to its projected growth in the next two decades 
(Richardson, 1984). These changes, coupled with the increased use of the Arctic as a 
deplovment area for the strategic forces of the Soviet Union, have forced the United 
States to develop an active and progressive Arctic research plan. The goals and objec- 
tives of this research are focused at solving scientific and technological problems related 
to the physical and biological components of the Arctic and the various processes that 
govern the behavior of these components. Increased knowledge in such areas as: the 
Arctic as a natural laboratory, the role of the Arctic in national defense and 
global/regional climate and weather are addressed in these objectives (National Science 
Foundation, 1987). 

The United States Navy’s interest in the Arctic has surged in the last several years. 
It has long been realized that the Arctic is an important naval operating area of the 
Soviet Union. The Northern Fleet is the largest of four Soviet fleets. It boasts 65% of 
the Soviet’s nuclear ballistic missile submarine (SSBN) force and 50% of its attack and 
cruise missile submarines (Polmar, 1986). Recent observation of Soviet fleet operating 
procedures and exercises seems to suggest that at an initial onset of hostilities a large 
majority of the Soviet naval forces would deploy relatively near the Soviet Union 
(Watkins, 1986). 

Advances in technology over the last ten vears have resulted in ranges up to 8300 
km for Soviet submarine-launched ballistic missiles. Development of submarine hull 
designs such as that of the Typhoon class, believed to be specifically designed for oper- 
ations under ice, gives the Soviets a significant strategic advantage. They now have the 
option of deploying their SSBN’s within Arctic waters and under the ice. This compli- 
cates the anti-submarine warfare (ASW) problem immensely and significantly reduces 
the vulnerability of these strategic platforms to detection by North Atlantic Treaty Or- 
ganization (NATO) ASW forces (U.S. Department of Defense, 1987). 
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Fig. 1. Arctic Ocean and adjacent seas: heavy dotted line is the defense early 
warning (DEW) Ine: large filled in circles are ballistic missile early warning 
sites (BMLEWS); heavy solid line shows the northern sea route; ice edge 
denotes average summer positions (Johnson, ev al., 1984). 


In the event of hostilities, it would also be in the interest of the United States to 
contain the Soviet fleet in the Arctic area. This requires a naval force able to effectively 
operate for extended periods of time in the Arctic arena at any time of the vear. Fig. 2 
shows the United States Navy's Arctic’/Cold Weather Surface Ship Plan developed by 
the Vice Chief of Naval Operations for Surface Warfare (OP-03). The current and 
projected severe weather operating capability for naval surface ships indicates the Navy's 
intention of expanding its operating limits to the verv edge of the marginal ice zone 
(MIZ) by the vear 2010. 
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Fig. 2. U.S. Navy’s Arctic/Cold Weather Surface Ship Plan. 


Our ability to operate in the seas bordering the Arctic Ocean is currently limited 
both by our knowledge of this area and by our abilitv to predict and avoid severe 


Weather. In response to the need for a better understanding of the Arctic environment, 


several major research efforts have been undertaken. The Marginal Ice Zone Exper- 
iment (MIZEX) series is one of these. 


B. THE MARGINAL ICE ZONE EXPERIMENTS 

The MIZ occurs geographically as a transition area between pack ice and the open 
ocean. It is an area of dynamic interaction involving the ice, the ocean and the atmos- 
phere. The MIZ boundary is dependent on the advection of ice from the north, 
advection of warmer water from the south and their interaction. Such interactions occur 
in the form of ocean waves and swell, mesoscale eddy grinding and from in situ ice for- 
mation and melt. The ice edge boundary varies seasonally by hundreds of kilometers. 
This covering and uncovering of vast amounts of water causes significant changes in the 
local atmospheric boundary layer (ABL) structure. In addition, the amount of open 
water has a direct bearing on the amount of moisture and heat available for exchange 
with the atmosphere. These fluxes of moisture and heat over large areas have a signif- 
icant impact in determining global climate (Johnson and Hawkins, 1987). 

The MIZEX series of experiments was developed to increase our understanding of 
the physical and biological processes of the MIZ. The measurement program was initi- 
ated with a small test phase in June and July of 1983. The following year a full-scale 
experiment was conducted from 18 May to 30 July 1984. This MIZEX involved the 
coordinated resources of eleven different nations including seven ships, eight remote 
sensing’meteorological aircraft and four helicopters. Data collected dunng MIZEX-84 
were used to study a variety of mesoscale physical processes involved in ice, ocean and 
atmospheric interactions (Johannessen and Horn, 1984). The operating area for 
MIZEX-84 is depicted in Fig. 3. 

A MIZEX-87 was conducted from 20 March to 10 April 1987 to collect data during 
the late winter regime. Fig. 4 shows the area of operations for MIZEX-87. The oper- 
ating area is in the same locale as that of MIZEX-84. The data obtained are significant 
in that this was the first major study of the Greenland Sea MIZ during the late winter 
period. Winter conditions are expected to produce substantially different dynamic bal- 
ances of momentum, heat and moisture as well as different interactions between ice- 
ocean conditions and their remotely sensed electromagnetic radiation properties 
(Johnson and Hawkins, 1987). One of the primary goals of the 1987 MIZEX field pro- 
grams was to provide the first data on important meteorological questions, including 
cyclogenesis and surface atmospheric boundary conditions in the winter MIZ (MIZEX 
Bulletin VIII, 1986). 
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Fig. 4. MIZEX-87 area of operations (Johnson and Hawkins, 1987). 





Three surface vessels were used to collect meteorological data during MIZEX-87. 
The Polar Circle, a Norwegian research vessel strengthened for movement in the ice, 
recorded three-hour surface observations, launched rawinsondes every six hours and 
made continuous measurements of mean wind, moisture, temperature, pressure and 
turbulent wind. The Polar Circle operated inside the West Fram Strait ice edge for the 


majority of the experiment between 74.5° to 79°N and 19°E to 6.5°W. Haakon Mosby, 


another Norwegian research vessel owned by the University of Bergen, Norway, oper- 
ated primarily at the ice edge or just off the ice edge and made measurements simular to 
those of the Polar Circle. The Ilaakon Mosby operated between 74.6° to 79°N and 
I9°E to 3.7°W. The Waldivia, a research ship fron: West Germany, recorded surface 
observations every three hours and launched rawinsondes every six ours. The Valdivia 
Operated in the open ocean away from the ice edge during most of MIZEX-87 between 
foto SO°Neand 17°E to 3°W. 

A MIZEX shore forecasting center was established at Tromso, Norway. At this 
center reports from the three ships were received and added to routine ship and shore 
station reports provided by the World Meteorological Organization (WMO) reporting 
network to provide analyses and forecasts during the experiment. Satellite imagery from 
the Advanced Very High Resolution Radiometer (AVHRR) mounted on the NOAA 9 
and 10 polar orbiters were used to help verify the analvses and track the observed 
weather svstems (Schultz, 1987). 


C. PURPOSE AND SCOPE OF THESIS 

As the operating tempos of militarv forces continue to increase in the Arctic, our 
knowledge of the environment and how 1 affects platform and weapons sensors will play 
a kev role in the deplowment of these assets. The ability to predict the occurrence of 
Various anomalies in electromagnetic (EM) propagation paths, such as radar ducts and 
radar holes or shadows, is a tactical necessitv. This thesis follows the work done by 
Willis (1987) who investigated the ABL of the summer MIZ and its effects on the 
refractivity of electromagnetic waves within this laver. Rawinsonde data collected by 
four ships during MIZEX-84 were used as input into the Integrated Refractive Effects 
Prediction Svstem (IREPS) to determine ducting parameters of height, thickness and 
strength. A spatial and temporal statistical analysis of the observed trapping lavers and 
associated ducts was completed. A simular analysis is performed in this thesis using 
rawinsonde data from the three vessels that participated in MIZEX-87. The late winter 
ABL 1s discussed and compared to the conditions observed during MIZEX-84. The 
spatial and temporal statistical analyses are compared to those in the MIZEX-84 study 
by Willis. A range-dependent ray trace model is used to see the effect of a horizontally 
varving refractive structure on the EM wave propagation pattern coming from a trans- 
mutter. Various refractive structures are used to simulate conditions observed during 
MIZEX-87. 


I. REFRACTIVITY IN THE ATMOSPHERE 


EM energy travels in the form of waves and can be classified according to wave- 
length or frequency. In the context of military applications the EM spectrum is often 
subdivided into categories of EM and electro-optical (EO) wavelengths. The EM por- 
tion of the spectrum includes wavelengths greater than I cm but less than 10 m. The 
EO portion of the spectrum is defined to be wavelengths from 0.4 to 100 um. This thesis 
directs its discussion only towards the EM portion of the spectrum. Further mention 
of EM radiation should be recognized to apply as such. 

The propagation of EM radiation through the atmosphere is controlled by the ver- 
tical gradients of temperature, moisture and pressure. The effect of these meteorological 
variables on the mean refractive structure of the atmosphere is most pronounced for 
wavelengths ranging from 1 cm to 10 min length. The frequencies associated with these 
wavelengths are classified in the radio spectrum as very high frequency (VHF; 30-300 
MHz), ultra high frequency (UHF; 300-3000 MHz) and super high frequency (SHF; 3-30 
GHz). In the radar spectrum these wavelengths equate to radar frequencies classified 
as the A-K radar bands. At longer wavelengths the ionosphere becomes the controlling 
influence. For wavelengths less than 1 cm but greater than 100 um atmospheric molec- 
ular absorption generally makes long-distance propagation impractical. 

EM radiation traveling in space at the approximate speed of 3 x 108m s~' experiences 
an effective resistance to its propagation speed as it passes through an atmosphere with 
varying constituents and properties. This resistance is caused by the impinging external 
electric field inducing polarization of the individual molecules in the medium. This 
polarization produces an opposing electrical field which has the effect of reducing the 
propagation speed. The opposing field is caused by orientation of the molecules with a 
dipole moment anti-parallel to the field or by distortion of the electron shell structure. 

The dielectric constant, €, is used to describe the interaction of EM radiation with 
a medium. The amount of interaction, and therefore the strength of the opposing field, 
is dependent both on the type of medium and on the frequency of the radiation. The 
relationship is given by (2.1), where c equals the velocity in a vacuum and v equals the 


velocity in the medium. 


Vv =c/ ell? (2.1) 


The index of refraction, n, is defined by (2.2) as the ratio of the velocity of an EM wave 
in a vacuum to its velocity in another medium and is equivalent to the square root of the 


dielectric constant. 


n=cv=el? ( 
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If an EM wave travels through a region of varving n it will be refracted towards the re- 
gion of higher n or towards slower wave speeds. This is a simple application of Snell's 


Law defined in (2.3) and illustrated in Fig. 5. 
n, sin @, =n, sin 0, (2.3) 


Generally, a denser material has a greater value of n. The density of the atmosphere 
usually decreases with height and thus EM waves are bent back towards the earth’s 
surface. 

The index of refraction in the atmosphere is a function of the value of pressure, 
temperature and water vapor pressure. The refractivity, N, defined in (2.4), is often used 


to represent refractive effects more conveniently. 
N =(n-1) x 10° (2.4) 


An equation to calculate the refractivity is given by Battan (1973) in (2.5) where P is 
pressure in millibars, T is temperature in Kelvin and e is water vapor pressure in milli- 
bars. This value of refractivity is simply a summation of the calculation of refractivity 


for dry air and for water vapor. 
N = 77.6P/T - 5.6e/T + 3.75 x 10°e/T? (2.5) 


Battan indicated that carbon dioxide also contributes to N but that its contribution is 
negligible at less than 0.1% of the total value. The refractivity of dry air is essentially 
the same over the full scale of the EM spectrum. With the addition of water vapor to 
the atmosphere, the refractivity develops a frequency dependence. This is due to the 
dipole nature of the water molecule and its differing response to varying frequencies of 
radio waves. The greatest effect is seen in the lower frequency EM waves identified 
earlier and results in a greater value of n than for higher frequency optical waves. 
When determining levels of EM transmission anomalies, the measure of the vertical 


gradient of refractivity is used. To account for the earth’s curvature, a modified measure 





Fig. 5. Graphical representation of Snell’s Law. 


of refractivity M, is defined bv Battan (1973) in (2.6) where N, is the value of N at a 


given height, Z. in kilometers. 
Mi Nea SZ (2.6) 


This modification aids one in the identification of trapping lavers and the associated 
ducts. 

A trapping laver is a vertical section of atmosphere where M decreases with height. 
EM ducting is a direct result of the presence of trapping lavers. An EM duct is a region 
in the atmosphere that tends to channel or confine EM waves. EM waves travelling in 
a duct incur lower spreading losses than in a normal atmosphere. This results in a signal 
strength that is greater at increased ranges and can lead to greatly increased radar or 
communications ranges. By knowing where levels of ducting exist, military tacticians 
can use such knowledge to plan aircraft attacks and surveillance routes. It can also be 
used to determine when counterdetection ranges of aircraft and surface radars are 
counterproductive to their useage. Finally, by knowing where enhanced and degraded 
ranges occur, one can better estimate initial detection ranges of hostile platforms trav- 
eling in or passing through these levels. 

Fig. 6 identifies three types of ducting situations. A surface based duct caused bv a 
surface trapping layer is depicted in profile (a). Profile (b) is an example of a surface 
based duct that occurs as a result of an elevated trapping layer. Finally profile (c) shows 
an example of an elevated duct caused by an elevated trapping laver. The thickness of 
the duct is determined by taking the height at the top of the trapping laver and sub- 


tracting from it the height at which an identical value of M is reached when drawing a 
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Fig. 6. Three types of ducting situations after Dotson (1987): profile (a) is a sur- 
face based duct caused by a surface trapping laver; profile (6) is a surface 
based duct caused by an elevated lJaver; profile (c) is an elevated duct 
caused by an elevated layer. 


vertical line down from the top of the trapping layer. This is best seen in Fig. 6c. Ifa 
transnutter is located within a duct the emitted waves will be concave upward with re- 
spect to the earth’s curvature as they travel through that portion of the duct with a 
positive M gradient. When the waves encounter the trapping laver they will be concave 
downward and curve back toward the earth. This ray theory concept is illustrated in 
Ric: 7. 
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Fig. 7. Possible ray paths in a ducting situation after Kerr (1951). 


The Navy's Integrated Refractive Effects Prediction System (IREPS) is a shipboard 
environmental data processing and display svstem used to predict the effects of refrac- 
uon on EM signals. The values of N and M gradients used by IREPS to define where 
trapping lavers occur are listed in Table 1 taken from Naval Ocean Systems Center 
(1981). When the vertical gradient of M is equal to 0 per km, an EM wave will be re- 
fracted downward with a curvature equaling that of the earth’s surface (Battan, 1973). 
If the M gradient has a negative value, the curvature will be greater than that or tne 
earth resulting in a ducting situation. The M gradient value of 0 per km 1s defined as 
the threshold value for trapping of an EM wave. 

Other anomalies of importance in atmospheric EM wave propagation are levels of 
super-refraction and subrefraction. The values of N and M gradients used by IREPS to 
determine levels of super-refraction and subrefraction are shown in Table 1. Super- 
refraction is defined to be where M gradients vary from 0 to 79 per km. In a super- 
refractive region the EM waves are bent towards the earth but not to the extent observed 
in a trapping layer. Radar propagation ranges are somewhat greater than normal in this 
region. In a subrefractive laver the EM waves bend away from the earth’s surface. The 
M gradient is greater than 157 per km. A subrefractive laver causes reduced radar ranges 
to occur. The values listed in Table 1 are also used in this study to determine levels of 
trapping, super-refraction and subrefraction. The ray geometry associated with the 


various conditions of refraction is depicted in Fig. 8. 
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Fig. 8. Ray geometry for various refraction conditions (Naval Ocean Systems Cen- 
ter, 1981). 
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Hl. THE MARINE ATMOSPHERIC BOUNDARY LAYER IN THE 
ARCTIC 


A. GENERAL 

A thorough understanding of the physical mechanisms influencing the depth of the 
atmospheric boundary layer (ABL) is essential in the study of atmospheric refractivity. 
Large mean gradients of temperature and moisture often occur at the boundaries of the 
ABL. As discussed in chapter two, these vertical gradients cause EM waves to refract 
according to the sign and strength of the refractivity gradient. This attenuation in the 
expected path can result in greatly increased or decreased wave propagation. 

The ABL is described by Dutton and Panofsky (1984) as a region of the atmosphere 
that varies in thickness from tens of meters in a stably stratified atmosphere to several 
thousand meters in an unstable case. All exchanges of momentum, moisture and ther- 
mal energy between the earth’s surface and the atmosphere take place within this layer. 
The ABL’s lower boundary is the surface of the earth or the ocean, and the flow within 
this laver is typically turbulent. 

Turbulence may be described as being either mechanically or convectively produced. 
The amount of turbulence observed in the ABL is a function of the surface roughness, 
the vertical gradient of the near-surface wind and the sensible and latent heat flux orig- 
inating from the ground and the atmosphere. Turbulence is produced mechanically from 
the flow of air over a rough surface by energy extracted from the gradient of the mean 
wind. A higher surface wind speed and a stronger vertical gradient of speed coupled with 
a rougher surface increases the amount of shear production of turbulent kinetic energy 
(TKE). This occurs to some degree regardless of the thermodynamic stability of the 
ABL. Sensible heating at the earth’s surface or the release of latent heat during cloud 
formation results in the convective production of turbulence. This causes the ABL to 
tend towards greater instability and increases both the level of available TKE and the 
mixing throughout the layer. As the amount of available TKE increases, the ABL is 
expected to increase in thickness (Wyngaard, 1973). 

The unstably stratified or convective ABL is subdivided into three distinct regions 
as depicted in Fig. 9. The lowest of these, the surface layer, ranges in height from a few 
meters to several tens of meters depending on the strength and direction of the thermal 


energy flux. Within this layer large vertical gradients of wind speed, temperature and 
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Fig. 9. The convective atmospheric boundary Jayer: mean vertical profiles of po- 
tential temperature (@), wind speed (UL), and specific humidity (Q) are 
superimposed. 


moisture occur. The strong vertical wind shear and the convection due to warmer sur- 
face temperatures produces the turbulence that transports moisture, heat and momen- 
hum Up into the nuxed laver. The nixed laver, directly above the surface laver, is 
relatively homogeneous in nature and composes the majority of the ABL. This laver 
approaches homogeneity as the turbulence produced within the boundary laver causes 
the erosion of local gradients that form within the mixed layer. The upper 50-100 meters 
of the convective ABL is depicted as a capping inversion layer. Within this laver the 
TKE available for mixing is consumed by a strong, increasing vertical temperature gra- 
dient. This gradient effectively prevents the exchange of physical properties with the free 
ecmosphere above (Wyneaard, 1973). 

A stable ABL is almost exclusively due to cooling at the surface and subsequent 
development of a positive vertical temperature gradient. Larger scale convectively gen- 
erated turbulence is suppressed in this situation, with very little mixing occurring above 
the surface laver. This results in a significant decrease in the exchange of momentum, 
heat and moisture between the earth’s surface and the atmosphere. 

The marine atmospheric boundarv laver (MABL) differs from the ABL over land in 
several ways. The ocean surface experiences a much smaller diurnal temperature vari- 
ation than that of the land. This causes the height of the MABL to vary relatively less 


diurnally. The ABL over land often experiences significant fluctuations in height due to 
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radiative heating and cooling. The humidity structure in the MABL is often much more 
pronounced than that over the land. The lower liquid boundary provides a ready source 
of moisture to the MABL. This results in geater vertical variations of moisture over the 
MABL surface layer and capping inversion (Dotson, 1987). The stability of the MABL 
is very nearly neutral over most of the world’s oceans. As air flows over the ocean for 
a long period of time its properties move toward equilibrium with the ocean surface. In 
this near-neutral state, variations in stability can have a significant effect on the surface 
fluxes of moisture, heat and momentum. Thus even small variations in the sea surface 


temperature can cause the MABL structure to change (Businger and Shaw, 1984). 


B. MARCH/APRIL SYNOPTIC CLIMATOLOGY AND OBSERVED 
CONDITIONS 
I. Unique Features 

Changing synoptic conditions can significantly affect the temperature, momen- 
tum and moisture distributions in the MABL and thus alter the level of mixing in and 
the height of the boundary laver. Several features that are unique to the Arctic region 
influence the synoptic weather patterns. These features are identified in Sater et al. 
(1971) as distinctive daylight and darkness regimes, the circumpolar vortex, surface snow 
and ice cover, and the flatness of the Arctic region. 

Distinctive regions of daylight and darkness are characteristic of the Arctic due 
to the tilt of the earth’s axis and the lower elevation of the sun. In the winter a perpetual 
darkness occurs and results in large radiational loss at the earth’s surface. During the 
summer, in conditions of nearly constant davlight, a net radiational loss is still present 
but not to the degree seen in the winter. These radiational losses cause cooling at the 
surface and development of a surface temperature inversion especially in the central 
Arctic. The months of March and April are a transition period between the two ex- 
tremes. The shift from total darkness to total daylight occurs over only a few weeks 
time. During MIZEX-87 the hours of davlight lengthened at the rate of 25 minutes per 
day from 9 hours of daylight in mid-March to 21 hours by mid-April (Schultz, 1987). 

The circumpolar vortex is the name associated with a large scale cold core low 
pressure region centered around the North Pole, and the belt of upper-level westerly 
winds that surround the low at its outer edge (Petterssen er al., 1956). Fig. 10 is a 
12-year mean of the 700 mb height field in decameters. The intensity of the circumpolar 
vortex is often associated with the packing of the 960 and 990 dm lines at this level. The 


surface cyclones and anticyclones that travel through the Greenland Sea region are 
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Fig. 10. Climatological 700mb heights in April from Sater et al., (1971). 


embedded in and steered by this upper level flow. Changes in the position and the 
strength of the vortex can radicaliy affect regional conditions (Sater ez al., 1971). The 
March, April climatology shows the vortex moving from its most southern position in 
the winter, when the packing and the strength of the vortex is the greatest, to its most 
northern position in the summer, when the vortex is weakest. 

The snow and ice cover so prevalent in the Arctic affects the meteorology of the 
area in several ways. The high surface albedo properties of snow and ice cause a much 
higher percentage of the incoming solar radiation to be reflected. This causes lower 
temperatures to be maintained over the ice compared to the ocean. The colder temper- 
atures reduce the moisture capacity of the air to such an extent that a large portion of 
the Arctic is often referred to as a desert (Sater ez al., 1971). The ice cover also effec- 
uvely blocks upward heat flux from the warmer ocean waters beneath it. This has a 


siguificant impact on the stability of the ABL over the ice. These areas tend towards 


greater stability and lower ABL thicknesses. The presence of an ice-ocean boundary 
creates areas of significant horizontal temperature gradients. These gradients have been 
closely correlated to instances of cyclogenesis, especially during conditions of off-ice 
wind flow (Petterssen er al., 1956). This factor is of great interest when considering the 
synoptic patterns in and around the MIZ. 

The flatness of the Arctic region promotes relatively unlimited flow of air into 
the region. Cyclones that are often generated along the polar front in midlatitudes are 
able to travel freely into the Arctic basin. Because of the low temperatures and lack of 
available moisture, cyclones frequently stagnate and fill in this region. Frequent cyclonic 
activity 1s encountered in the Greenland Sea MIZ as cyclones from the midlatitudes are 
steered by the circumpolar vortex up into this region. 

2. Spring Climatology 

The region in which MIZEX-87 was conducted is under the influence of several 
different climatological surface pressure features. Fig. 11 shows the 12-year 
climatological mean sea-level pressures for April. The synoptic pattern around the 
MIZEX area is characterized by a weak and variable surface high over northern 
Greenland and a belt of low pressure extending from the southern tip of Greenland, 
through Iceland and up through the Norwegian Sea to the northern tip of Norway. The 
belt of low pressure corresponds to cyclones moving through this region that are 
strongly influenced by upper-level steering. The wind fields in the MIZEX area are 
generally easterly due to the presence of the migratory lows to the south (Sater er al., 
1971). 

The climatological movement of cyclones is more clearly shown in Fig. 12. The 
primary and secondary cyclone tracks indicate that, in April, cyclones move through the 
Greenland Sea area in a northeastward direction. The primary storm track originates 
off the east coast of Iceland and moves up to a region off the northern coast of Norway. 
The secondary track follows the southeast coast of Greenland northward and then cuts 
across the Greenland Sea, passing just south of the Svalbard Islands and into the same 
sink region as the primary track. 

3. Meteorological Conditions Observed in MIZEX-87 

The synoptic and mesoscale meteorological features observed during MIZEX-87 
have been examined and separated into five distinct periods by Schultz (1987). During 
period one, from 20-23 March, the MIZEX area was under the influence of a high 


pressure system centered over the Fram Strait. Period two was from 24-27 March. 
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Fig. 11. Climatological sea-level pressures for April from Sater et al., (1971). 


During this time the area was influenced by a mesoscale boundary layer front that 
formed just west of the Svalbard Islands. The Fram Strait high weakened during this 
period and surface pressures fell. The third period, from 28 March to 31 March was 
influenced by a small scale low that persisted until the last two days and then began to 
fill. The fourth period was from 1-3 April. During this period the Fram Strait region 
was under the influence of a small high pressure system. A synoptic scale low moving 
northeast passed to the south of the MIZEX area but did not significantly influence the 
local conditions. The final period. from 4 to 10 April, saw the area under the influence 
of two separate synoptic-scale lows. These lows formed about 50 km off the east coast 
of Greenland at about 71°N and 15°W. The first low moved across the Greenland Sea, 
following the secondary storm track depicted in Fig. 10 and stalled about 100 km 
southwest of the Svalbard Islands on 5 April. The second low caught up and merged 


with the first low on § April creating a larger, more intense system that significantly 





Fig. 12. Climatological tracks of cyclones in April from Sater et al., (1971): Pri- 
marv tracks indicated by solid lines. Secondary tracks indicated by dashed 
lines. 


affected the boundary laver structure of the region. These five periods will be examined 
in greater detail in chapter five to determine how the observed meteorological features 
altered the refractive structure of the atmosphere. 

The observed synoptic and mesoscale features were consistent with the 
climatology of the MIZEX area discussed earlier. The region was dominated bv high 
pressure except during the movement of synoptic lows along the secondary cyclone 
track. The occurrence of the boundarv-laver front is also a common feature of this re- 
gion. Winds tended to be generally from the north during MIZEX-87. This differs from 
the climatological easterly winds proposed by Sater er al. (1971) suggesting a lessened 


influence by the transiting low pressure systems on the wind direction during this period. 


C. MIZEX-87 BOUNDARY LAYER STRUCTURE 

The structure of the ABL during MIZEX-87 varied considerably both spatially and 
temporally. Variations in the height of the ABL have been correlated by Guest e¢ al. 
(1988) with respect to the synoptic situation, wind speed and direction relative to the ice 
edge, and horizontal location with respect to the ice edge. In the absence of turbulence 
measurements, the level of mixing in the boundary layer was determined by visual in- 
spection of the rawinsonde data. Levels where the potential temperature was nearly 
constant with height were considered to be well-mixed. This is strictly true only when 
no liquid water drops are present. The potential dewpoint temperature and specific hu- 
midity were used to identify cloud layers and well-mixed regions respectively. Two hu- 
midity variables were used in determining cloud lavers and mixed layers because no 
single humidity variable can easily show their location (Guest and Davidson, 1988). 

Investigation of the rawinsonde data revealed that the boundary laver normally 
presented a well-mixed appearance with a pronounced capping inversion when the area 
was under the influence of high pressure. Fig. 13a shows a plot of potential temper- 
ature, potential dewpoint temperature and specific humidity for a rawinsonde launch 
from the Valdivia taken at 0608 UTC on 23 March. During this time the Valdivia was 
located just south of a high centered over the Fram Strait. The boundary laver is well- 
mixed up to 650 m, where a capping inversion separates it from the free atmosphere 
above. During periods of lower pressure and frontal activity as in periods 2 and 5, the 
boundary laver increased in thickness to heights which often exceeded 1500 m. Fig. 13b 
is an example of this situation taken from the Haakon Mosby at 1147 UTC on 6 April. 
The layer is well mixed up to 1600 m and the inversion is nonexistent or very weak at 
best. Lindsay er al. (1986) presented similar results in demonstrating that the inversion 
tends to descend in response to warming conditions in the free atmosphere associated 
with higher pressure and rises during periods of cooling and greater instability brought 
on by low pressure systems. 

Correlation of the observed surface wind speed with the boundary laver height is 
depicted in Fig. 14. The median height of the top of the mixed layer generally increases 
as wind speed increases in a somewhat linear fashion from a median value of 500 m at 
1.5 ms-! to a median value of 1075 m at 13.5 ms“. As discussed earlier, increases in 
the surface wind speed aid in the generation of TKE available for mixing of the bound- 
ary layer. However, higher wind speeds were experienced during periods of lower pres- 


sure, and this almost certainly causes some contamination of the correlation. 
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Fig. 13. Vertical profiles from two different pressure regimes: potential temper- 
ature (bold solid line); potential dew-point temperature (dashed line); 
specific humidity (solid line). 


Wind direction and horizontal position with respect to the ice edge also appears to 
show some correlation with the height of the MABL. Fig. 15 from Guest ez al. (1988) 
shows that on-ice winds produced deeper well mixed layers than off-ice winds. The on- 
ice Winds however. had much more individual scatter in the layer heights than did the 
off-ice winds. This scatter is not shown in Fig. 15. Guest reported that the on-ice winds 
often produced multiple inversion lavers as the offshore air flowing onto the colder ice 
surface produced new boundary lavers. These new boundary layers were associated with 
lower inversion heights, however the median change 1s not nearly as drastic as that ob- 
served in the off-ice flow. 

The off-ice flow shows median mixed laver heights to be below 200 m over the ice 
and then to increase rapidly as the air moves out into the open ocean. This increase can 
be attributed to the increase in heat flux from the ice to the open ocean resulting in a less 
stable more turbulent boundary laver over the ocean. Fairall and Markson (1987) de- 
termined that the heat flux near the ice edge was close to zero and that the heat flux over 
the ocean was 30 W m=. The increased heat flux over the ocean results in an increase 


in the energy available for mixing in the layer and entrainment of the free atmosphere. 
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Fig. 14. Variation of median mixed layer height with surface wind speed. 


This entrainment will cause the inversion height to increase while decreasing the inver- 
sion strength. 

Winds parallel to the ice edge with the ice to the right of the wind flow showed a 
rather mixed pattern in the median mixed layer heights. Guest ev al. (1988) suggests that 
the zig-zag type pattern may be a result of varying surface roughnesses normal to the 
wind vector or across the ice edge. This varying roughness can induce vertical motions 
due to frictional convergence or divergence and thus cause variations in the level of 


nuxing. 


D. COMPARISON OF MIZEX-84 AND MIZEX-87 MABL STRUCTURE 

Several differences were noted in the boundary layer structure between the 
MIZEX-84 summer case and the late winter, early spring MIZEX-87 data. Fig. 16 from 
Guest er al. (1988) shows the median Well mixed layer heights for both MIZEX-84 and 


MIZEX-87 as a function of horizontal position. Both data sets displayed a similar 
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Fig. 15. Miedian mixed layer heights from Guest et al., (1988): off-ice, on-ice and 
parallel winds are indicated by triangles, squares and circles, respectively. 


pattern, however the mixed layer depth was consistently higher for the springtime case. 
It was also noted that the strength of the inversion laver was typically stronger during 
the summer MIZEN. In addition, the air above the inversions was warmer over 
ice-covered areas aS opposed to open ocean areas. This suggested a greater degree of 
subsidence over the ice. This variation in upper air temperatures was not evident during 
MIZEX-87, and the boundary layer appeared to be much less stable. Overland (1985) 
deternuned that the winter boundary laver in the MIZ was generally shallower with 
sharper inversions. The summer regime was marked by a deeper boundary layer with 
weaker capping inversions. The data presented by Guest er al. (1988) do not seem to 
follow this trend. There are several possible reasons for the height differences observed 
during the two MIZEX’s. One possible reason is that the cyclone tracks during 
MIZEX-87 were farther north than those in MIZEX-84, which may have reduced the 
divergence in the area. Another cause suggested is that during the summer the warm 


air associated with the nearby land masses may have been advected into the atmosphere 
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Fig. 16. Comparison of mixed layer heights from Guest et a/., (1988): MIZEX-84 
is indicated by squares, MIZEX-87 is indicated by triangles. 


over the inversions above the ice by an ice breeze circulation causing stronger inversions. 
This supply of air is not present over the cooler continents during the March’April time 
frame (Guest ez al., 1988). Another possible cause in the height differences may be due 
to more periods of onshore flow during MIZEX-87 as compared to MIZEX-84. 
Onshore flow resulted in higher median mixed laver heights, and this could explain part 
of the difference. One other possibility is that the colder air temperatures experienced 
during MIZEX-87 would have created greater surface-to-air heat fluxes in the 
MIZEX-87 boundary layer. Greater heat fluxes would influence the turbulence driving 


the nuxing in the layer and would yield higher inversion lavers. 
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IV. DATA ACQUISITION AND PROCESSING 


A. DATA ACQUISITION AND SOURCES OF ERROR 

The data used to analyze the refractive conditions observed during MIZEX-87 were 
obtained from the three surface ships operating in and around the MIZ and from the 
Norwegian Meteorological Institute based in Tromso, Norway. The research vessels 
Polar Circle, Haakon Mosby and Valdivia launched rawinsondes at approximately 
6-hour intervals and collected meteorological data using ship-mounted instruments. 
Visual surface observations were also made on each ship every three hours. 

1. Polar Circle 

The Polar Circle collected data from 22 March through 11 April 1987. The ship 

was equipped with a meteorological instrument station mounted on a platform located 
16 m above sea level and just forward of the low mast. The station failed to give accu- 
rate measurements for wind speeds above 9 ms-'!, During these episodes wind speed 
was measured by a sonic anemometer located next to the meteorological station. For 
wind speeds less than 9 ms-! the two instruments showed excellent agreement with 
maximum speed differences of only 0.3 ms-!. The instrument accuracy of the winds at 
these slower speeds should be less than or equal to this difference. Distortion of the 
winds by the ship’s structure were estimated to be less than 5% due to the instruments 
being located away from any blocking structures. An exception to this is for winds 
coming directly from astern when errors may have been as great as 20%. Measurements 
of wind direction were subject to errors in the magnetic compass of the meteorological 
station due to magnetic anomalies caused by the ship or by the local geology. This re- 
sulted in extreme directional errors in the winds for short periods of time. For a majority 
of the time the wind direction measurements from the meteorological station on Polar 
Circle were estimated to be accurate to within 20 degrees. The surface temperatures 
were estimated to be accurate to within 0.2°C and the relative humidity to within 5°%% 
of the stated value. The above data were sampled continuously onboard the Polar Circle 
at a sampling rate of 12 seconds. Surface atmospheric measurements presented in this 
thesis represent 10-minute averages of these data. Pressure was recorded using auto- 
matic pressure sensors. These sensors are not as accurate as other instruments but 
provided increased data coverage. Errors in pressure readings were not greater than 2 
mb (Guest and Davidson, 1988). 


2. Haakon Mosby 
The Haakon Mosby measured the same parameters as the Polar Circle using the 
same instrumentation as the Polar Circle except for wind speed. A miniature cup 
anemometer was used on the Haakon Mosby for higher wind speeds. Errors in wind 
direction Were not as extreme because the ship did not cause significant magnetic dis- 
tortion. Wind directions were usually accurate to within 5°. Winds from the stern of 
the ship exhibited errors of up to 20% in speed and direction due to flow distortion. 
The other instrument accuracies are the same as those reported for the Polar Circle 
(Guest and Davidson, 1988). 
3. Valdivia 
The Valdivia measured the same parameters as the previous two ships, however, 
due to timing errors in the recording instruments the most reliable data were taken dur- 
ing the surface observations every three hours. The accuracies of the measurements 
Were similar to those reported for the Haakon Mosby (Guest and Davidson, 1988). 
4. Rawinsondes 
Rawinsondes manufactured by the VIZ corporation were used on the Polar 
Circle and Valdivia to obtain vertical profiles of atmospheric parameters. The 
VIZ-Beukers W-8000 RP+ Microsonde system and its associated software were used 
on these two ships to obtain measurements of pressure, temperature, dew-point tem- 
perature, relative humidity, altitude, wind speed and wind direction. The sampling in- 
terval of the rawinsonde was every 10 seconds for a period of 50 minutes. This provided 
a maximum of 300 vertical data points per launch. The wind speed and wind direction 
Were measured by tracking a signal from the rawinsonde using an Omega tracking sys- 
tem. The temperatures measured using this system were accurate to 0.2°C and the 
dew-point temperature was accurate to 1.0°C. Heights were calculated to be accurate 
to within 30 meters. The wind speeds were accurate to | ms-~! except in regions of 
strong wind shear where the accuracy dropped to 2 m s-! (Guest and Davidson, 1988). 
The Haakon Mosby used a rawinsonde system manufactured by the Vaisala 
Corporation. Omega tracking was used for calculation of wind speed and direction. The 
measured parameters were identical to those measured by the VIZ system. 
A total of 169 successful launches were made from the three ships. The Polar 
Circle recorded 66 launches from 24 March to I1 April, 1987. Valdivia recorded 48 total 
launches from 19 March to 3 April, 1987. Haakon Mosby recorded 55 launches from 
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24 March to 11 April, 1987. Of these, two of the launches resulted in only partial data 
sets. 

The vertical gradient of moisture, as noted previously, is a key variable in de- 
termining the refractivity gradient dM'dZ. Willis (1987) described a possible error in the 
dewpoint curves of 37% of the MIZEX 84 rawinsonde profiles examined in her study. 
In the suspect profiles the dewpoint temperature equaled the air temperature just below 
the inversion, indicating saturation and the presence of a cloud layer. At the inversion 
base the dewpoint normally decreases with height. This decrease is due to the drving 
out of the air caused by subsidence in the air above the inversion. The profiles examined 
by Willis did not follow this pattern. The dewpoint curve was observed to remain satu- 
rated through the inversion and decreased sharply at a variable height above the inver- 
sion. Fig. 17 from Willis (1987) shows an example of this phenomenon. 

Willis offered several physical explanations as to the cause of the suspect pro- 
files. One suggestion was the possibility that warmer and more moist air, originating 
over the northern boundary of the Gulf Stream, could have been lifted and advected into 
the Arctic region. This would require the winds to be primarilv from a southerlv direc- 
tion which was not the case during MIZEX 84. Another explanation proposed by 
Ohtake et al. (1982) attributed the saturation above the inversion layer to the increased 
water vapor flux from open leads in the ice during sub-zero temperatures. This would 
cause the atmosphere to be supersaturated with respect to the ice. This proposal as- 
sumes that the measurements are being taken in an area predominantly covered by ice. 

A second possibility addressed by Willis was that of a humidity sensor problem. 
It was suggested that the humidity sensor may have been lagging after encountering an 
area of saturation. This could have been caused by moisture accumulating on the sensor 
in the form of condensation or possibly frost/ice. It was noted that the dewpoint often 
decreased rapidly shortly after ascending through 0°C into warmer temperatures. It is 
possible that frost formed on the humidity sensors would be able to melt and subse- 
quently evaporate in the drier air above the inversion. The effect of such an error in the 
moisture gradient would be to increase the height of the base of a trapping layer as well 
as lessen the strength and the thickness of the layer. Willis applied a correction to the 
suspect profiles to eliminate the saturated dewpoint curve along and above the inversion. 
Her results were reported using both the corrected and uncorrected data sets. The fre- 
quency of ducting events using the corrected data was found to be only half that meas- 


ured using the original data. 
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Fig. 17. Suspect Dewpoint Curve from Willis (1987). 
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The dewpoint problem described by Willis was observed very infrequently in the 
MIZEX-87 rawinsonde profiles. In 13 of 169 launches (7.7%) a situation similar to that 
described by Willis in the MIZEX-84 data was noted. Of the 23 incidences of ducting 
observed, five of them were affected bv this phenomenon. Of the 13 suspect profiles, five 
were onlv near saturation as thev followed the temperature curve up through the inver- 
sion. In 10 of the 13 profiles the rapid decrease in the dewpoint occurred before the 
temperature reached 0°C. Additionally, 10 of the 13 profiles were launched from the 
Polar Circle and occurred while the Polar Circle was inside the ice edge or very near to 
it. This fact lends support to the theory presented by Ohtake er al. (1982). Although 
an error in the humidity sensor could be a problem, it was determined that the dewpoint 
problem was not extensive enough in the MIZEX-87 data to warrant the correction. 

Future experiments in the MIZ should address this issue specifically to deter- 
mine if the phenomenon is real or due to instrument error. Willis (1987) suggested that 
rawinsonde profiles could be compared to the dewpoint profiles of an aircraft descending 
from the dry, warm air above down through the inversion. This would provide more 


conclusive evidence as to the authenticity of the strange dewpoint trace. 


B. DATA PROCESSING 

Rawinsonde data serve as the primary data source for this study. Intitial processing 
of the data included checking the raw data for any errors caused by loss of signal, loss 
of tracking, multiple readings at the same level or balloon bursts. These errors were 
identified automatically by the system software. A second review of each data set was 
then done to manually identify any errors missed by the software. An interpolation 
routine was used to smooth the data through the error points. In certain cases wind 
speeds and directions were not available due to a loss of the rawinsonde tracking signal 
by the Omega tracking system. This did not affect refractivity measurements. Surface 
data collected from shipboard instruments was plotted in time series of temperature, 
relative humidity, pressure, wind speed and direction. 

A Fortran computer routine was prepared to calculate the refractivity N and the 
adjusted refractivitv value M defined in chapter two. Equations 2.5 and 2.6 were used 
to calculate these values. The vapor pressure required in equation 2.5 was calculated 
using the formula in (4.1) after Buck (1981), where e is vapor pressure in mb, P is pres- 


sure in mb and TD is the dew-point temperature in K. 


e = (1.0007 + (3.46 x 107° x P)) x 6.1121 x expo? 9027 PI2%0.974TD)) (4.1) 
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Values of the vertical gradients of M were determined between each level. These values 
were compared to the categories listed in Table | to determine if the atmosphere between 
two consecutive levels was either subrefractive, normal, super-refractive or trapping in 
nature. The gradient values were subsequently used to determine the strength and the 
thickness of individual EM ducts. The height of a duct was defined as the opumum 
coupling height (OCH). This is the height at which dM ‘dZ changes from a positive to 
a negative value. The height of super-refractive and subrefractive lavers was defined as 
the height at the base of each layer. 

A range-dependent ravtrace program, developed by Wavne Patterson of the Naval 
Ocean Systems Command (NOSC), was used to see the effects of a horizontally varving 
refractive structure on the raypaths of EM waves (Patterson, 1987). Refractivity profiles 
similar to those observed in the previous case studies were input into the program at 
range intervals of 50 km. The height of a source transmitter was varied to see how the 
raypaths would be altered when the source was placed within a ducting layer as opposed 
to outside of one. The transmitter was given a beam width of 1.5° and an elevation 
angle of 0° with the horizontal. The program incorporates linear ray theory into the 


determination of the ray paths and is frequency independent. 
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V. TEMPORAL STUDIES 


Incidences of EM ducting, super-refraction and subrefraction were studied to deter- 
mine variations in the frequency of occurrence and changes in individual heights of the 
lavers as a function of time. A synoptic regime case study covering the period 20 March 
through 10 April, 1987 was performed to see how the refractive structure was influenced 
by the synoptic-scale forcing during each period. A diurnal study was performed to see 


if the height and frequency of refractive events varied over a 24 hour period. 


A. SYNOPTIC REGIME STUDY 

A synoptic regime study was performed by separating the MIZEX-87 data set into 
a series of five periods based on the synoptic situation present during each period. These 
periods were defined and described in depth by Schultz (1987). 

1. Description of the Synoptic Periods 

Period one covered the first four days from 20-23 March, 1987. Fig. 18 shows 
the surface pressure pattern for 0000 UTC 20 March, 1987. A surface high of 1045 mb 
is situated over northeast Greenland. <A surface low of 990 mb is present between 
Iceland and southern Norway. During the four-day period the surface high moved 
northeastward and then southward over the Fram Strait and began to weaken. The low 
pressure to the south moved towards Iceland and began to weaken. Fig. 19 shows a high 
of 1035 mb situated over the Fram Strait and a 1000 mb low just to the east of Iceland 
at 1500 UTC 23 March. The prevailing winds were easterly during this period. The 
Valdivia recorded wind speeds that varied from a maximum of 24 m s~! on 20 March to 
no wind on 23 March. During this time the Valdivia was moving westward from the 
Svalbard Islands towards the MIZEX-87 operating area. 

The second period was from 24-27 March, 1987. During this four-day period a 
boundary layer front developed about 100 km west of the Svalbard Islands. The 
boundary layer front is a mesoscale phenomenon observed to form along the ice edge 
boundary when cold dry air is advected from over the pack ice to the open water. The 
air over the water warms due to sensible heat flux from the water. Cloud formation en- 
hances this warming with the addition of latent heat. The front develops at the bound- 
ary between this modified air and the unchanged air mass over the ice or land surface. 
The front moved rapidly westward across the Fram Strait on 25 March and stalled upon 
reaching the MIZ. Sea-level pressure in the Fram Strait dropped from 1035 mb to 1017 
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Fig. 18. Sea-level pressure analysis, OUU0 UTC 20 March 1987. 


mb over the period. An upper level high continued to persist in association with the 
front. As the front passed the three vessels, wind speeds increased by about 10 m s-! and 
wind direction shifted clockwise from northerly to southeasterly. Fig. 20 shows that the 
surface low pressure area associated with the boundary laver front is centered just south 
of the Fram Strait on 1200 UTC 25 March. 

During the third period, from 28-31 March, 1987, development of a weak sur- 
face low was observed over the Fram Strait. The sea-level pressure pattern for 1200 
UTC 28 March is depicted in Fig. 21. The pressure in the MIZEX-87 operating area 
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to a minimum of about 1010 mb on 29 March and then rose steadily throughout the 
period to 1024 mb on 31 March. During this period the upper level high continued to 
weaken. Winds were primarily from the north and northeast during the period and 
ranged from 3 to 10 ms-'. By 3! March the high pressure to the north showmin Fig. 
22 began to have the dominating influence on the Fram Strait region. 

Period four was from 1-3 April, 1987. During this period a synoptic-scale low 
formed just to the east of Iceland. Fig. 23 shows an inverted trough analyzed over the 


Fram Strait region with high pressure to the east and west. Although the low was fairly 
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large, its effects were not significantly felt by the ships operating in the Fram Strait re- 
gion. The pressure in the MIZEX-87 operating area rose from 1024 mb on | April to 
1032 mb on 3 April as the high pressure on either side of the strait dominated the 
large-scale forcing. Winds during the fourth period were generally from the north. A 
significant shift in the wind direction occurred at 0400 UTC on 3 April from northerly 
to southeasterly. The wind speed was constant at about 7 ms~! before the shift and 
became highly variable after the shift, ranging from 3 to 13 ms-'. This shift in wind 


direction Was probably due to the passing of a weak front associated with the low 
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Fig, 21. Sea-level pressure analysis, 1200 UTC 28 March 1987. 


pressure center located to the south as depicted in Fig. 24. The Valdivia departed the 
MIZEX-87 operating area on 2 April en route to Tromso, Norway. 

The fifth period included the last seven days from 4-10 April, 1987. Two more 
major cyclones formed about 400 km north of Iceland. On 5 April the first low moved 
to the northeast from a position 50 km off the east coast of Greenland to a position 150 
km south of the Svalbard Islands. The final position of this low is shown in Fig. 25. 
The second low formed in the same area as the first low and followed its track, over- 


taking and merging with it on 8 April at a position 200 km southwest of the Svalbard 


36 





Fig. 22. Sea-level pressure analysis, 0000 UTC 31 March 1987. 


Islands. The pressure in the MIZEX-87 operating area remained fairly constant around 
1015 mb through § April. The winds were generally from the north although shifts in 
wind direction to the southeast were measured on 6 Apnil along with an increase in wind 
speed of 5 ms-'. The Polar Circle and Haakon Mosbv departed the MIZ late in the 
evening on 8 April and on 9 April encountered a low pressure system formed from the 
merger of the two previous lows. The Haakon Mosby recorded a drop in pressure from 
1015 mb to 1003 mb over 12 hours time. Wind speeds increased from 7 to 14m s~' and 


a shift in wind direction from the north to the southeast was noted. After moving 
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through the low, the observation on Haakon Mosby showed a pressure increase back 


to 1015 mb and a shift in wind direction back to the north. The wind speed decreased 
to about 6-8 m s-!. 


Refractive Structure During Each Period 
a. Period One (20-23 March) 


During period one the Valdivia was just arriving on station in the MIZ area. 
The Polar Circle and Haakon Mosby were still en route to the area. On 20 March, the 


Valdivia Was moving northward towards the Svalbard Islands. Fig 18 shows the position 
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of Valdivia at 0000 UTC. By 23 March the Valdivia had moved to the west, nearer to 
the MIZ. During this period the Valdivia recorded increasing surface pressure as it 
moved towards the high pressure centered over the Fram Strait. Fig. 26 shows that 
during this period one duct was detected by the Valdivia. This duct was detected at a 
height of 709 m and was caused by the boundary laver inversion forming at that level. 
Analvsis of the rawinsonde profiles in Fig. 27 reveals that the inversion continued to 
strengthen and the mixed laver shallowed as the ship moved into an area of increased 


subsidence aloft. 
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Fig, 25. Sea-level pressure analysis, 1200 UTC 5 April 1987. 


Fig. 28 shows the super-refractive laver occurrences and heights over time. 
Super-refractive layers were observed in seven of the nine rawinsondes launched during 
this period. In five of the seven cases the super-refractive layer was caused by the pri- 
miarv inversion. The other incidences of super-refraction occurred due to sharp drops in 
the dewpoint temperature with height as the rawinsonde passed into a drier laver not 
associated with the capping inversion . 

Subrefractive layers occur primarily above the boundary layer in areas 


where there is either a substantial negative vertical temperature gradient or, more 
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importantly, a large positive moisture gradient. This often occurs where differential 
advection of moisture into an area above the boundary laver produces the necessary 
moisture gradient. If the moisture is sufficient to form a cloud, the subrefractive laver 
will be just below cloud base. Fig. 29 shows that, in period one, subrefractive layers were 
observed in three of the launches. On 21 March the subrefractive laver was associated 
with a thin moist layer at about 800 m. Multiple subrefractive lavers were detected in 
the other two launches. Of these, one was detected just below a cloud layer at 115 m. 
Another one of the multiple subrefractive layers was based at the surface and had a 
thickness of 75 m. I[t was associated with a relatively drv surface relative humidity of 
42% that increased over 75 mto 65%. The large variation in relative humidity may be 
due to an inaccurate value for the first data point measured. The first few data points 
of a rawinsonde profile are often observed to be in error and strong low-level gradients 
such as those above may be suspect. The final three subrefractive layers were caused 
by thin layers of more moist air above the well-mixed layer. 
b. Period Two (24-27 March) 

The refractive structure of the second period was affected significantly by 
the inversion that had developed over the MIZ during the previous period. Ducts were 
observed during eight of the 16 launch episodes. A launch episode includes all of the 
rawinsondes launched by the three ships at approximately six-hour intervals. All but 
one of the ducting incidents occurred during the first two days of the period. Fig. 26 
shows the duct heights relative to time. In each case the duct was caused by the pres- 
ence of a strong capping inversion over the well-mixed layer. Fig. 30 shows that the 
capping inversion remained relatively strong around the MIZ through 26 March. At 
2300 UTC on 26 March the boundary laver front described in the synoptic overview of 
period two passed over the Polar Circle. The rawinsonde profile for 1125 UTC 27 March 
in Fig. 30 indicates that the front has destroved the inversion and caused the temperature 
to be fairly well-mixed up past 2000 m. Similar trends were observed in the profiles of 
the other two ships as the boundarv laver front advanced westward across the Fram 
Strait. Only one duct was recorded during the last two days of the period at an elevation 
of 1229 m where the inversion was beginning to reform following frontal passage. 

Incidents of super-refraction closely paralleled those of ducting during the 
second period. A super-refractive layer was observed by at least one of the ships during 
13 of the 16 launch episodes. Of the 27 super-refractive layers identified, five of them 


Were adjacent to a trapping laver. 17 of the super-refractive layers were caused by the 
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Fig. 27. Rawinsonde profiles from synoptic period one. 
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Super-refractive layer heights relative to day and time of launch . 
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Subrefractive layer heights relative to day and time of Jaunch. 


mig. 29, 


boundary layer inversion. Fig. 28 shows the scatter of the super-refractive layers during 
period two. Similar to the trapping lavers, most super-refractive layers were detected 
before the boundary layer front moved through. The super-refractive layers observed 
at the Valdivia and Haakon Mosby about 24 hours after frontal passage indicated that 
the inversion layer was beginning to reform but was not strong enough to produce con- 
ditions required for the formation of trapping layers. 

Very few subrefractive lavers were seen during period two. Over the four- 
day period, Fig. 29 indicates that only four episodes produced subrefractive lavers. The 
scarcity of subrefraction during this period is due to the high levels of moisture 
throughout the boundary layer. Very few significant increases in the moisture gradient 
were observed in the profiles. The passage of the boundary layer front increased the 
mixing to higher levels. Although there was little indication of upper level advection of 
moist lavers into the region, any gradients that were forming would have been reduced 
by the enhanced mixing. 

c. Period Three (28-31 March) 

During period three the MIZEX-87 region was under the influence of a 
weak low pressure system. The increased instability associated with this low and a fur- 
ther weakening of the upper-level high pressure resulted in only one ducting episode at 
the very end of the period as shown in Fig. 26. This duct occurred as the pressure had 
steadily risen from a low of 1010 mb on 29 March to 1024 mb by 31 March. Fig. 31 
shows that the boundary layer was rather deep throughout the period and lacked a 
well-defined capping inversion. The inversion is just begining to develop in the last 
profile shown in Fig. 31 at a height of 1250 m. 

Numerous detections of super-refractive layers were noted in period three 
during 12 of the 16 launch episodes. Of the 25 layers observed in the rawinsonde profiles 
of the three ships, 18 were associated with a weak inversion above the boundary layer 
accompanied by decreases in the dewpoint. The weaker inversions were caused by a 
reduction in the subsidence rate due to the weakening upper-level high pressure. The 
remaining super-refractive layers were detected just above upper-level bands of increased 
moisture where the moisture gradient dropped significantly. These super-refractive lay- 
ers did not occur in coincidence with a strong temperature inversion. 

Subrefractive events were noted during eight of the 16 launch episodes. The 
subrefractive layers formed at levels where a transition to a substantially more moist 


layer was observed. In only three of the events was saturation apparent, indicating 
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Fig. 30. Rawinsonde profiles from synoptic period two. 
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cloud formation. This suggests that although cloud layers are often observed just above 
subrefractive layers, the presence of clouds is not a mecessary requirement for 
subrefration to occur. 

d. Period Four (1-3 April) 

The fourth synoptic period exhibited a similar refractive structure to that 
observed during the end of period one and the beginning of period two. High pressure 
dominated the MIZEX area and resulted in the development of a stronger capping in- 
version above the well-mixed layer. The stronger inversion provided the necessary M- 
gradient for ducting to occur during seven of the first eight launch episodes. Fig. 32 
shows that the duct heights varied between 750 and 1100 m as the inversion height 
fluctuated over the period. Four rawinsonde profiles in Fig. 33, from the Haakon 
Mosby, indicate that the inversion remained pronounced until early afternoon on 2 
April. By late afternoon on the same day the inversion began to weaken and it appears 
that a deep layer of mid-level moisture was advected into the area. The wind shift noted 
in the synoptic discussion of period four signaled the passage of the front associated with 
the low to the south. The effects of the advected moisture and the frontal passage are 
evident in the final profile of Fig. 33 which shows the boundary layer nearly saturated 
with no evidence of an inversion. 

Super-refractive events observed during the fourth period followed a pattern 
similar to that observed in the ducting episodes. Super-refractive events were recorded 
in all of the launch episodes during the first two days of the period. No events were re- 
corded during the last day. The super-refractive layers identified in Fig. 34 were associ- 
ated with a temperature inversion in 16 of the 37 reported layers. Nine more of the 
super-refractive layers occurred just above layers of clouds. The remaining layers oc- 
curred where the dewpoint decreased sharply over a very short distance through a series 
of narrow bands of increased moisture. 

Subrefractive layers were detected primarily during the first two days of the 
period as indicated in Fig. 35. A majority of these layers were found at levels where the 
moisture gradient was strongly positive but did not reach a point of saturation. These 
subrefractive layers were found between the thin super-refractive layers mentioned 
above. 

e. Period Five (4-10 April) 
Period five was the longest (7 days) of the synoptic regimes. Refractive 


events of interest occurred during 4-5 April and 7-8 April. By 4 April the sea-level 
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Fig. 31. Ravwinsonde profiles from synoptic period three. 
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Fig. 35. Subrefractive layer height relative to day and time of launch. 
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Fig. 36. Rawinsonde profiles from synoptic period five. 


54 


pressure had become steady at 1015 mb following the passage of the front on the day 
before. Surface temperatures dropped almost 10°C on 4 April as the winds shifted to the 
west. northwest bringing in colder air from the Greenland land mass. The atmosphere 
tended towards more stable conditions and the boundary layer inversion began to re- 
form. Just outside the MIZ a rawinsonde profile showed two inversions in the vertical. 
The lower inversion formed at about 300 m. The second inversion capped a well-mixed 
layer at about 1200 m. Fig. 36 shows the development of the inversions during 4 April. 
The inversions were not strong enough to produce more than one ducting event as in- 
dicated in Fig. 32, however, a significant number of super-refractive lavers were created 
by the two inversions. These are plotted in Fig. 34. The super-refractive lavers found 
between 500-1000 m appeared to occur just above the tops of some deep cloud lavers 
that descended nearly to the surface. Incidents of subrefraction were scarce during 4-5 
April due to a lack of moist layers in the profiles. 

During 6 April the MIZEX area experienced relatively few refractive events. 
This was probably due to increased convergence and thus mixing in the MABL brought 
on by the second low pressure system of the period as it moved across the Greenland 
Sea. Numerous refractive events were detected during 7-8 April. Fig. 32 shows that the 
Polar Circle recorded three low-level ducts and two higher ducts over the ice while the 
Haakon Mosby, located outside the ice boundary, recorded three ducts at about 650 m. 
Fig. 34 shows a large number of super-refractive lavers detected on these two days dur- 
ing the first six launch episodes. The layers were found at nearly the same levels as the 
ducting. Subrefractive events during these two days were located just below cloud layers 
both above and below the inversion levels. 

3. Trends in the Refractive Structure. 

Several trends were noted as a result of the regime study. Refractive events 
occurred with much greater frequency during periods of increased atmospheric stability. 
During the periods dominated by high pressure the boundary layer capping inversion 
was able to develop and strengthen, providing the conditions necessary for trapping 
layers to occur. During frontal passages or when the MIZ was under the influence of 
low pressure, the vertical stability of the atmosphere was degraded. This results in in- 
creased vertical mixing and reduced subsidence in the air above the MABL and a sub- 
sequent decav of the inversion. 

Many more super-refractive layers than trapping layers were observed during 


MIZEX-87. In addition, the bases of the super-refractive layers were higher than the 
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duct heights. The average duct height recorded during MIZEX-87 was 835 m. The av- 
erage super-refractive layer height was 1048 m. The difference in these heights was noted 
during every series of launch episodes in which ducts and super-refractive layers were 
observed. This pattern is due to the way the inversion develops. Upper-level subsidence 
causes adiabatic warming and drving of the air as it descends. As it comes in contact 
with the well-mixed layer a temperature inversion and a negative moisture gradient begin 
to form. As the strength of these gradients increase, the M-gradient required for super- 
refraction occurs before that for trapping. As subsidence increases, the rate of descent 
of the air above the MABL exceeds the upward growth of the layer due to entrainment 
generated by mechanical and convective turbulence. This results in a depression of the 
inversion to a lower level and a strengthening of the inversion gradient. If the 
subsidence is strong enough a trapping laver will eventually develop at this lower level. 

Another observation of interest was that in all cases a ducting event was asso- 
ciated with both a temperature inversion and a negative moisture gradient. Both of 
these conditions had to be present at the same level before a trapping layer would form. 
Super-refractive lavers were frequently observed when only one of the two conditions 
was met. A majority of the time that condition was a sharp decrease in moisture. 

Subrefractive layers were almost always seen when individual cloud layers were 
present above the boundary layer. However, saturation was not a requirement for sub- 
refraction. Many instances of subrefraction were recorded when clouds were not pres- 
ent. During periods of increased instability the number of subrefractive lavers detected 
dropped significantly due to increased mixing of the moisture profile. Subrefractive 
events were attributed primarily to sharp positive increases in the moisture gradient 
rather than negative vertical temperature gradients. 

4. Comparison to MIZEX-84 Regime Studies 

Ducting episodes were studied and described by Willis (1987) over the period 
19 June - 14 July, 1984. Six different synoptic regimes were determined based on the 
large-scale synoptic forcing that influenced the area. Fig. 37 shows the percentage of 
launch episodes during each regime in which a ducting event was observed. Both cor- 
rected and uncorrected data results are given. For the uncorrected case, regimes three, 
four and six recorded ducting in greater than 50° of the launch episodes. These regimes 
were identified by Willis as being under the influence of ridging caused by high pressure 
systems located over Greenland and Norway. The subsidence aloft produced by this 


ridging caused the inversion to develop sufficiently during these regimes to provide 
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conditions necessary for persistent ducting to occur. Regimes one, two and five were 
dominated by low pressure systems moving through the Greenland Sea and resulted in 
a reduced number of ducting episodes. These results are consistent with the ducting 
observed during MIZEX- 87. 

In the analysis of EM refractivity around the Greenland Sea MIZ, it is apparent 
that the large-scale forcing has a dominant effect on the upper boundary of the MABL 
and on the strength of the associated inversion. A majoritv of the ducts and super- 
refractive events observed during MIZEX-87 were associated with this inversion. The 
large-scale synoptic forcing must be considered as a primarv factor in the development 


of conditions necessary for anomalous refractive events to occur. 


B. DIURNAL STUDY 
A study was conducted to determine if the diurnal cycle had any effect on either the 
heights of refractive events or on their frequency of occurrence. The diurnal cycle was 
driven by the large change in the number of daylight hours over the three week period. 
The number of daylight hours increased from 9 to 21 hours over the course of 
MIZEX-87. Figs. 38, 39 and 40 are plots of the heights of ducting, super-refractive and 
subrefractive events respectively versus the time of day the event was detected. 
1. Diurnal Variations in Height 
Duct heights and super-refractive layer heights would be expected to follow 
diurnal height changes in the MABL inversion. A thorough study of the rawinsonde 
profiles during each diurnal cvcle was completed. No discernable diurnal patterns in the 
MABL height were obvious. Dhiurnal changes would be expected due to radiational 
heating and cooling at the cloud tops below the base of the inversion. Diurnal! changes 
in the fluxes responsible for generating the mixing in the boundary layer could also cause 
variations in the boundary layer height. Time series of the sea-surface temperature, 
surface air temperature, relative humidity and surface winds were examined. Although 
large-scale fluctuations caused by the synoptic forcing were apparent, diurnal fluctu- 
ations Were not observed in these time series. Rawinsonde profiles were examined over 
24 hour periods to see if the changes in the heights followed any pattern. Diurnal vani- 
ations were not observed. Movement of svnoptic-scale systems through the area, as well 
as continuous movement of the ships could be causing any diurnal variations in the 
refractive layer heights to be effectively masked. The scatter observed in the heights of 
all three diurnal plots was significant during each launch period and no appreciable 


trends Were evident when comparing one launch period to another. 
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2. Diurnal Variations in Frequency of Occurrence 

Examination of the diurnal scatter plots seemed to indicate a larger occurrence 
of refractive events during the midday and early evening launch periods. A statistical 
study was completed to determine if these variations in frequency of occurrence were 
significant. The results of this study are presented in Table 2. The number of 
rawinsondes launched during each launch period varied. In order to account for this in 
the statistical study, a percentage of occurrence for each of the refractive effects was 
calculated for each period based on the total number of rawinsondes launched during 
that period. 

Incidents of ducting were detected more frequently during the midday and early 
evening launch periods compared to those for early morning or late evening. This ob- 
servation is significant in that the trend is opposite to that normally observed or ex- 
pected. Normally conditions for ducting tend to improve during the hours of darkness 
due to radiational cooling at the top of clouds. This causes an increase in the temper- 
ature jump across the inversion which should lead to a greater occurrence of ducting at 
night. An explanation for the increased frequency of ducting during the midday and . 


early evening hours is not readily apparent. This observation could be of significant 


tactical interest if it is observed in other data sets and should be investigated further. 










TABLE 2 
PERCENTAGE OF OCGLRRENCE OF REFRACTIVE Ey Ena 
BASED ON TOTAL LAUNCHES PER EAUNCH CYCEE 


Poenius ot 1) Rega 
Lavers Lavers 
0500:0600 UTE 






Super-refractive lavers followed a pattern simular to that of ducting although the 
highest percentage of super-refractive layers was 64°%o, observed at midday. It may be 
possible that the increase of ducting rates between the nudday and early evening periods 
is related to the decrease in super-refractive rates during the same times. if the inversion 
strengthened during the interval between the two periods, ducting conditions could exist 
where previously super-refractive conditions existed. This suggestion, while plausible, 
was not readily apparent in the data set. Investigation of this observation was compli- 
cated by frequent changes in the large-scale forcing present over the region and by con- 
tinuous movement of the ships. Further studies of this phenomenon should be made 
with data measured from a fixed position during a period of relatively stable large-scale 
forcing. 

Subrefractive events showed little diurnal change in frequency of occurrence. 
A midday percentage of 40% was higher than the rate of the other three periods which 


Were nearly equal at 30%. The cause of this is not obvious. 
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VI. SPATIAL STUDIES 


Several spatial studies were conducted to investigate variations in EM refractivity 
over horizontal distances. In the first study, variations in EM ducting, super-refraction 
and subrefraction were examined based on position relative to the ice edge. The second 
study looks at cases where rawinsondes were launched at approximately the same time 
to determine how homogeneous the refractive structure of the atmosphere 1s over hori- 
zontal ranges. A range-dependent ray tracing program was used to show how inhomo- 


geneitv in the refractive structure affects EM wave propagation. 


A. HORIZONTAL VARIABILITY RELATIVE TO THE ICE EDGE 

Considerable variability in the heights of the various refractive layers was observed 
relative to the ice edge. Variations in duct thickness and strength were also noted. The 
results of these findings are summarized in Table 3. The ice edge boundarv was taken 
to be where 50% of the surface was covered bv pack ice. This boundary remained fairly 
constant during MIZEX-87. The average position of the ice edge is indicated by the 
dashed line in Fig. 4. For the purpose of this study the area within 25 km of either side 
of the ice boundary was considered to be within the MIZ. This boundary is for statis- 
tical purposes only as the actual width of the MIZ varies considerably depending on the 
season, currents and wind direction ‘speed. 

Guest and Davidson (1988) determined that the well-mixed layer height was greater 
over the ocean than over the pack ice. This was described in chapter three and illus- 
trated in Fig. 16. In chapter five the duct heights and super-refractive laver heights were 
associated with the capping inversion at the top of the well-mixed layer during many of 
the launch episodes. One would expect the duct and super-refractive layer heights to 
follow a pattern similar to that of the well-mixed layer relative to the ice edge. This 
trend was observed during the statistical study to a certain extent. To more closely ex- 
amine the variations in height, thickness and strength, the data were plotted in scatter 
plots as a function of distance from the ice edge. Fig. 41 is a plot of duct height versus 
distance from the ice edge. The mean and median values of duct height were calculated 
over the ice, MIZ and water. The duct heights increased from a mean value of 632 m 
over the ice to 890 m over the water. The median values for these two areas showed a 
similar trend. The duct heights over the MIZ were the highest with mean, median values 
of 982 m and 994 m respectively. The increase in the heights over the MIZ is an 
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unexpected feature and was also observed for super-refractive and subrefractive layers. 


Several explanations for this observation will be discussed later in this chapter. 


TABLE 3 
VAR Bibi iy IN REFRACTIVE STRUCTURE 
RELATIVE 10 IME ICE EDGE 
(values given as mean; median in meters) 


eee Over Ice “larginal Over Water 
Variable Ice Zone 
Duct , 
Height 632/530 982/994 890/793 


Duct x ihe | 
Duct . . | 


Super-refractive pee . | 
Laver Height 620/753 1332/1217 1192/1061 


Subrefractive ie - | 
Laver Height 624/728 1644/1782 1369) 1511 


* Duct strength measured in M-units 





Duct thickness is plotted relative to the ice edge in Fig. 42. The thickness of the 
ducts was relatively constant over the ice and varied less than 10 m for 80% of the ducts. 
Over the ice the average duct thickness was 42 m. In the MIZ the thickness varied 
greatly although the mean and median values closely resembled those measured over the 
ice. The duct thickness over the water was highly variable, ranging from less than 20 
m to greater than 100 m. The mean and median values of 54 m and 59 m were somewhat 
larger than those over the ice or MIZ. 

Duct strength was higher over the water than over the ice or the MIZ. Fig. 43 in- 
dicates that over the ice no ducts were detected with a strength greater than 3 M-units. 
Over the MIZ the ducts were even weaker. Over the water the strength of the ducts 
varied from values that approached zero to duct strengths as high as 7 M-units. 70% 
of the ducts detected over the ocean had strengths less than 3 M-units. 

Super-refractive layer heights and subrefractive layer heights were also examined 
relative to the ice edge. Fig. 44 is a scatter plot of super-refractive laver heights versus 


distance from the ice edge. Super-refractive laver heights over the ice resembled those 
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observed for ducting. The mean super-refractive laver height over ice was 620 m which 
is slightly less than observed for ducting. The median super-refractive layer height of 
753 m over the ice is greater than the median value of 530 m for duct heights and more 
clearly represents the observed situation. Over the MIZ, the super-refractive layer me- 
dian height nearly doubles to 1217 m and then falls slightly over the water to a median 
height of 1061 m. Subrefractive layer heights, illustrated in Fig. 45, followed the same 
general pattern observed for super-refractive layers. A median height of 728 m over the 
ice increased to 1782 m over the MIZ and then decreased to a median value of 1511 m 
over the water. 

To more graphically show the tendencies of the heights relative to the ice edge two 
plots-were constucted. Fig. 46 shows the first plot. In this plot all incidences of ducting 
are shown when at least two ships launched rawinsondes within the same launch epi- 
sode. If a duct was observed by more than one of the ships in a given launch episode 
the lines between the ducts were connected. Some evidence of an increasing slope going 
from right to left is apparent in this plot. The second plot, shown in Fig. 47 is similar 
to the first. However, in this plot super-refractive layers are added to the trapping layers. 
The super-refractive layers had to be associated with the same meteorological condition 
that produced the duct during that episode. In most cases this was the primary inver- 
sion. Fig. 47 indicates the positive slope even more clearly and supports the results ob- 
served using the mean and median values. 

The spatial variation in the heights of the various refractive layers relative to the ice 
edge is partly a function of the strength of the surface fluxes available for the formation 
of TKE in the MABL. The wind flow during most of MIZEX-87 was northerly. This 
sets up a flow pattern that brings relatively cold, dry air off from the pack ice over the 
warmer ocean surface. Air that remains over the pack ice is not influenced by increased 
surface heat and moisture fluxes and maintains a shallower mixed layer. This is con- 
sistent with observations by the Polar Circle as it operated within the pack ice. Air 
flowing over the warmer ocean surface experiences an increase in both heat and moisture 
fluxes. This causes increased turbulence and results in a deeper well-mixed layer. The 
higher refractive events detected over the water were often associated with the strong 
temperature and moisture gradients found at the top of the well-mixed layer. 

Over the MIZ higher mean and median heights were observed for all of the refrac- 
tion categories discussed. Based on the above physical discussion, one would expect 


these heights to fall between those observed over the ice and water. There are several 
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Duct height relative to distance from the ice edge. 
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Duct thickness relative to distance from the ice edge. 
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Fig. 44. 
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Fig. 46. Slope of duct heights recorded during the same launch episode relative to 
distance from the ice edge. 
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possible explanations for the deviation of the heights over the MIZ. Mean and median 
values of the heights included all detections of ducting and super-refraction and not only 
those associated with the inversion over the well-mixed layer. Multiple ducts in the 
vertical were recorded five times during MIZEX-87. Multiple super-refractive layers 
were observed during 28 different launches. In these cases only the lowest of the layers 
was associated with the primary inversion. A recalculation of the mean values of height 
was done using only those refraction events that could be associated with the inversion 
layer of the MABL. The mean values for duct heights showed no significant change over 
the ice. Over the MIZ the mean height dropped from 982 m to 901 m. Over the water 
the mean dropped from 890 m to 772 m. This does not contribute to an explanation of 
the greater height values over the MIZ. A more probable explanation for the higher 
duct heights is that the rawinsonde launches within the MIZ region were conducted 
during times when the MABL was relatively deeper than what it was over the total 
MIZEX-87 period. Further investigation confirmed this suggestion. Of the seven ducts 
detected within the MIZ area, five of them were detected at times when the boundary 
laver inversion was relatively weak and the height of the mixed layer was higher than 
normal. The recalculation of the super-refractive layer mean heights resulted in a drop 
in the mean height over the MIZ from 1332 mto 916 m. Over the water the mean 
dropped slightly from 1192 mto 1117 m. The recalculated mean heights reflect the in- 
crease in the well-mixed laver heights observed by Guest and Davidson (1988). 

Subrefractive layers were most closely correlated to areas just under cloud bases or 
under bands of higher moisture. Multiple levels of subrefraction were not observed over 
the ice. Conversely, over the MIZ and water there were 26 launch episodes that re- 
corded multiple levels of subrefraction in the vertical. These layers ranged from the 
surface up to over 3000 m in height. Most of these episodes occurred in open water 
within 50 km of the ice edge. 

The results of this study show that during the late winter/early spring the presence 
of ice has a definite effect not only on the height of the MABL but also on the level at 
which ducting, super-refraction and subrefraction are observed. Over the ice refraction 
anomalies were observed at relatively lower heights and were singular events. Over the 
MIZ and water the heights increased and greater incidences of multiple refractive layers 
were observed. 

Willis (1987) performed a simular statistical analysis of EM ducting relative to the ice 


edge. A comparison of results indicates that ducts were observed much more frequently 
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in the summer during MIZEX-84. Willis reported that ducts were detected during 
48.4°o of the launches in the uncorrected data and 30% in the corrected data. This is 
compared to a detection rate of 21% during MIZEX-87. Although the frequency of 
ducting events differed substantially, the mean values of duct height were reasonably 
close to the uncorrected MIZEX-84 values. Willis reported an uncorrected over-ice 
mean height of 502 m and a height of 848 m over the water. These compare to the re- 
spective values of 632 m and 890 m observed during MIZEX-87. After the correction 
for dewpoint error was applied, the mean duct heights over the ice and water were low- 
ered to 370 m and 560 m for the MIZEX-84 data. The correction applied to the 
dewpoint curves in the MIZEX-84 study caused the gradient of moisture near the in- 
version to be weaker. It also lowered the height at which the dewpoint began to de- 
crease, to one that more closely matched the base of the inversion. Based on the method 
in which the correction was applied and on the results of the MIZEX-87 study, it is 
thought that the frequency of occurrence of ducting 1s more accurate in the case of the 
uncorrected data for MIZEX-84. Conversely, the average duct heights relative to the ice 
edge are best represented by the values determined from the corrected data. 

The lower mean heights observed during MIZEX-84 could be due to the tendency 
towards stronger and more frequent ridging over the Fram Strait during the summer 
season caused by intensification of the Greenland High. This would result in a stronger 
subsidence rate aloft and a more depressed well-mixed layer height. The greater air/sea 
temperature differences observed during MIZEX-87 would also have increased the sur- 
face heat flux during the spring season. This would cause a relative increase in the TKE 
available for mixing in the boundary layer and a subsequent increase in the inversion 
base. 

Average duct thicknesses were two to three times greater during MIZEX-84 than 
during MIZEX-87. A similar trend was observed in the average duct strengths. These 
differences would be expected due to the stronger inversions that develop during the 
summer. The inversion strength is weaker in the spring due to lower subsidence rates 
compared to the summer. This effect is partly offset by the greater surface heat flux 
expected during the spring season. The heat flux causes increased mixing in the bound- 
ary laver and increased entrainment of the air above the inversion. This not only raises 
the base but also increases the temperature jump across the inversion. The large-scale 
forcing, causing the rate of subsidence to vary, appears to have a much more significant 


effect on the strength and thickness of the ducts. 
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B. VARIABILITY OF THE REFRACTIVE STRUCTURE BETWEEN SHIPS 

Of the 169 successful rawinsonde launches, 130 met the requirement for this portion 
of the spatial study of being launched within one hour of at least one other ship’s 
rawinsonde. These rawinsondes were launched during 53 seperate launch episodes. The 
spatial variability between ships was examined for each type of EM refractivity. Table 
4 gives a summary of the results. 

1. Statistical Results 

In 17 of the 53 launch episodes at least one duct was detected. Of these 17 ep- 

isodes, one resulted in ducts being reported by all three ships. A duct was reported by 
two of the ships in five of the launch episodes. The remaining 11 launch episodes in- 
volved a duct being reported by only one ship. The specific circumstances surrounding 
several of these launches will be discussed in three case studies presented later. Of the 
26 individual ducts detected, 66% were associated with the boundary layer capping in- 


version. 


TABLE 4 
PERCENTAGE OF DISTRIBUTION OF REFRACTIVE 
LAYER MULTIPLICITY PER LAUNCH PERIOD 


29% % 





Super-refractive events were noted in 44 of the 53 launch episodes. All three 
ships reported super-refractive lavers during five of the launch episodes. Two ships re- 
ported super-refractive layers during 21 of the episodes and only one ship reported a 
super-refractive layer during the remaining 18 episodes. 46% of the super-refractive 
layers detected were associated with the primary boundary layer inversion. In the case 
studies super-refractive layers detected on one ship are related to ducts detected on an- 
other ship during the same launch episode. 

Subrefractive layers were noted during 38 of the 53 launch episodes. Of these, 
three ships reported concurrent subrefractive layers during only two of the episodes. In 


nine of the episodes subrefraction was reported by two of the ships. The remaining 27 
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episodes had only one ship detect a subrefractive layer. The large number of single ship 
episodes seems to suggest that subrefractive events mav occur over smaller horizontal 
areas and may be caused by smaller-scale phenomena such as individual cloud layers or 
narrow bands of increased moisture levels. 
2. Case Studies 

A series of case studies was completed on various ducting episodes in order to 
more clearly understand the cause of the inhomogeneity observed in the refractive 
structure of the atmosphere during MIZEX-87. 

a. Case One 

Rawinsondes were launched from all three ships at approximately 1130 
UTC 24 March. Six hours later the Valdivia and the Haakon Mosby launched again. 
Polar Circle did not launch at the 1730 UTC launch period. Fig. 48 shows the positions 
of the three ships at 1200 UTC and 1800 UTC with respect to each other and the ice 
edge. In this and the following figures showing ship’s position, the Polar Circle is iden- 
tified by a circle, the Haakon Mosby by a square and the Valdivia by a triangle. At both 
launches the Valdivia was about 350 km to the north of Haakon Mosby with the Polar 
[ircle between. 

At the 1130 UTC launch the Haakon Mosbv and Polar Circle were less than 
20 km apart and in the open water. The area was under the influence of the Greenland 
High and winds were off-ice. All three ships detected ducts during the first launch. Fig. 
49 shows the potential temperature, dewpoint and specific humidity curves for the three 
rawinsondes. The associated M and N profiles for these rawinsondes are shown in Fig. 
50. The thicker line represents the M-profile and the thin line the N-profile. The ducts 
were detected at 579 m on the Valdivia, 968 m on the Haakon Mosby and 710 m on the 
Polar Circle. The strange dewpoint profile on the Haakon Mosby could be in error. A 
duct located at the base of the temperature inversion would have a height of 830 m. 
Even with this correction there is still a significant difference in the duct heights at a 
range of less than 20 km. 

At the 1730 UTC launch the temperature and dewpoint profiles in Fig. 51 
indicate a weaker inversion that has been lifted bv about 100 m over the past six hours 
at the Valdivia’s position. This is reflected in the duct height as it increased to 678 m. 
This weakening was probably caused by a slight warming of the well-mixed layer and a 
reduction in subsidence aloft as the Valdivia approached the boundary layer front de- 


veloping to the west. Conversely, the inversion detected on the Haakon Mosby has 
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Fig. 48. Ship positions at 1200 UTC and 1800 UTC 24 March 1987. 


decreased in height and strengthened. The duct height decreased to 818 m and the duct 
thickness increased from 42 to 96m. It1is apparent that the Haakon Mosby ts still under 
the influence of the high pressure svstem during the 17350 UTC launch. 

In this case study the large scale features of a well-mixed layer and a 
capping inversion existed throughout the MIZEX area but there was considerable vari- 
ability in the heights of the ducts detected by each ship. 


b. Case Two 
In this situation a condition was observed in which the Polar Circle, posi- 


tioned in the ice, recorded a duct at 1005 m. The Valdivia, stationed about 500 km to 
the southeast, recorded a super-refractive laver at 892 m. The Haakon Mosby was po- 
sitioned midway between the other ships and showed no evidence of ducting or super- 
refraction. Fig. 52 shows the positions of the three ships at 1800 UTC 26 March as 
well as the potential temperature, dewpoint and specific humidity profiles for each ship 


within one hour of that time. The M and N profiles for this case are given in Fig. 53. 
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Fig. 49. Rawinsonde profiles for launch episode of 1130 UTC 24 March 1987. 
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M and N profiles for'launch episode of 1130 UTC 24 March 1987. 
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Fig. 51. Rawinsonde protiles and the associated MI and N profiles for the launch 
episode of 1730 UTC 24 March 1987. 
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A closer look at the synoptic picture during this launch episode explained 
the strange refractive structure observed during this case. The Haakon Mosby observed 
frontal passage of the boundarv laver front at about 1600 UTC 26 March. The front 
effectively dissipated any inversion that had been present before and normalized the 
refractive structure. The Valdivia had experienced the frontal passage about 36 hours 
earlier and a redevelopment of the inversion is evident in Fig. 52 although it has not 
strengthened enough to permit ducting. The Polar Circle was still under the influence 
of the stable conditions aloft above the MABL. Observation of the launches just prior 
to and after the 1708 UTC launch by the Polar Circle revealed that the inversion base 
was lifting and the inversion was Weakening 1n response to the approaching boundary 
laver front. 

This case shows how a mesoscale meteorological phenomenon such as a 
boundary laver front can significantly alter the refractive structure of the atmosphere 
over a relatively short distance. 

c. Case Three 

This case examines four consecutive launch periods from 1130 UTC I April 
to 0530 UTC 2 April. Fig. 54 shows the relative positions of the ships during this period. 
The Valdivia was positioned over 400 km to the south of the other two ships. During 
this period the Valdivia was in an area under the influence of a passing synoptic-scale 
low and no incidents of ducting or super-refraction Were detected. At the same time the 
Polar Circle, inside the ice edge, was alwavs within 50-100 km of the Haakon Mosby, 
located just off the MIZ in the open water. Both of these ships were under the influence 
of high pressure. 

During the 1130 UTC launch on | April, the Haakon Mosby detected a 
duct at 871 m. The Polar Circle, stationed 100 km to the west and in the ice, recorded 
a super-refractive layer at 851 m. At the 1730 UTC launch, the duct detected bv the 
Haakon Mosby had descended to 822 m and the super-refractive layer observed six 
hours earlier by the Polar Circle had developed into a duct at 869 m. By the 2330 UTC 
launch time the ships had moved to within 60 km of each other. The Haakon Mosby 
detected a duct at 901 m and the Polar Circle observed a duct at 801 m. On the final 
launch at about 0530 UTC 2 April, the duct heights had increased on both vessels. 
Haakon Mosby detected a duct at 1017 m and the Polar Circle recorded a duct at 898 
m. It appears that the inversion associated with the ducts has lifted at about the same 


rate between the third and fourth launch episodes. 
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Fig. 52. Ship positions at 1800 UT 26 March 1987 and rawinsonde profiles for the 
1730 UT launch episode. 
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Ship positions from 1200 UT 1 April to 0600 UT 2 April 1987. 


Fig. 54. 


Figs. 55-58 show the rawinsonde profiles and their associated M and N 
profiles for the eight rawinsonde launches in which refractive events were observed. It 
can be seen in these figures that the duct heights increased as the inversion laver began 
to lift during the last two launch episodes. This was probably a result of the approach- 
ing front of the low pressure system to the south causing an increase in the instability 
in the boundary layer and a decrease in the subsidence aloft. The slope in the duct 
heights was positive going from ice to water during three of the four launches. 

In this case it was observed that for ships relatively close to each other ( less 
than 100 km ) and under the influence of the same large-scale forcing, the vertical 
movement of the ducting layers was fairly consistent. The individual heights of the ducts 
were observed to vary from 50-100 m at each launch episode. Over the horizontal ranges 
in question this amount of variation appears to be quite small. The effect of these var- 
lations on EM wave propagation will be discussed shortly. 

3. Comparison to MIZEX-84 Results 

The refractive structure of the MABL examined by Willis (1987) using 
MIZEX-84 data also indicated a lack of horizontal homogeneity. Ducts were detected 
by at least two ships during the same launch episode only 20 of 61 times for the uncor- 
rected data. This number dropped to only four after the dewpoint correction had been 
applied. Willis attributed some of the variations in duct height during individual epi- 
sodes to the dewpoint saturation problem. However, the point was also made that in- 
homogeneity in the refractive structure seemed to be the norm rather the exception in 
the data. 


C. EFFECTS OF HORIZONTAL INHOMOGENEITY IN THE REFRACTIVE 
STRUCTURE 

Variability in the horizontal refractive structure is of significant importance to any- 
one interested in exploiting that structure for tactical purposes. It has been shown that 
the refractive structure around the MIZ varies greatly over range and is highly depend- 
ent on both the large-scale forcing of the atmosphere and the more isolated effects of the 
ice cover and the ice/ocean interface. IREPS is the only operational system currently 
used bv the U.S. Navy to analyze the effects of the refractive structure on radar propa- 
gation. The input into IREPS is based on one sounding and the assumption is made 
that the atmosphere is horizontally homogeneous about 85% of the time (NOSC TD 
659, 1981). Dotson (1987) demonstrated the necessity to obtain high resolution data to 


effectively identifv small-scale refractive variations in the vertical. He also provided 
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Fig. 55. Rawinsonde profiles and associated Mand N profiles for 1130 UTC launch 
episode on | April 1987. 
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Fig. 57. Rawinsonde profiles and associated NI and N profiles for 2330 UTC launch 
episode on | April 1987. 
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Fig. 58. Rawinsonde profiles and associated M and N profiles for 0530 UTC launch 
episode on 2 April 1987. 
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evidence that horizontal inhomogeneity in the atmosphere is often more pronounced 
between two platforms than was previously assumed. This inhomogeneity must be taken 
into account when making predictions using IREPS which assumes a homogeneous 
refractive structure. A firm understanding of the forces driving the atmospheric condi- 
tions in between any two positions of interest 1s essential to the useful implementation 
of refractive effects prediction systemis. 

A range-dependent ray tracing program developed by Patterson (1987) was used to 
describe the effects of a varying refractive structure on the propagation paths of EM 
waves. Five different refractive structures were selected as input into the model. These 
included a horizontally homogeneous refractive structure, a refractive structure in which 
the trapping laver was allowed to slope upward with range, a refractive structure that 
began with a well defined trapping layer and changed to a normal refractive structure, 
a normal refractive structure that developed into a trapping layer and a trapping layer 
that abruptly changed into a normal refractive structure at a short range. Each of these 
refractive structures will be examined separately and compared to each other. 

1. Refractive Structure One 

A horizontally homogeneous refractive structure was used for the first study. 
The refractive structure was normal with the exception of one trapping laver based at 
S00 m. This height was chosen as it is approximately the average height of all the ducts 
observed during MIZEX-87. The thickness of the trapping layer was 50 m which was 
also close to the average thickness. 

Fig. 59 shows the ray path for a transmitter positioned at the base of the trap- 
ping layer at $00 m. An increased concentration of EM wave energy is identified where 
the trace is shaded darker. The effect of the duct created by the trapping layer is evident 
as the wave energy is focused within the duct creating a dark narrow bundle of rays 
propagating straight out from the transmitter. Another effect of the trapping layer is to 
generate a radar shadow zone or radar hole just above the layer. This is a zone which 
is void or nearly void of EM wave energv due to the trapping laver’s refracting the en- 
ergy downward. Radar holes are of great tactical importance because a contact flying 
into such an area would become effectively invisible to the tracking svstem. It is noted 
in Fig. 59 that the radar hole begins to form at a range of about 40 km. 

Two more runs were made with the homogeneous refractive structure but with 
the transmitter located above the trapping layer. The first of these depicted in Fig. 60 


has the transmitter at an altitude of 900 m. Several changes are readily apparent in the 


on 


ray trace after the height change. The narrow beam of focused rays directed straight 
ahead is no longer present. With the transmitter above the trapping layer the wave 
fronts intersect the trapping layer at an angle of incidence that exceeds the critical angle 
necessary for the duct to trap the energy. The radar hole is more distinct and the hole 
begins to form at about 30 km. Additional runs demonstrated that the radar hole would 
form closer to the transmitter as the transmitter height was reduced to a level no less 
than the top of the trapping laver. As the transmitter was allowed to increase in height, 
the radar hole formed further away from the transmitter. Fig. 61 shows the radar hole 
forming at 75 km with the transmitter at a height of 1200 m. When the transmitter was 
raised to a height of 1600 m the lower rays of the beam pattern were trapped in the duct. 

In the final run of this refractive structure, the transmitter was placed below the 
trapping layer at 700 m. Fig. 62 illustrates this situation. A radar hole did not form in 
this case. An area of increased energy is evident just above the trapping layer and slopes 
slightly upward. 

2. Refractive Structure Two 

This situation represents one in which a trapping layer, associated with the 
MABL capping inversion, is increasing in height as the range increases. The slope of the 
height increase was 100 m every 50 km. The trapping layer thickness was maintained 
at 50 m. This situation is similar to that observed during MIZEX-87 and represents the 
slope observed between observations over the ice and over the water. The initial 
refractivity profile in this series is identical to the one used in the previous homogeneous 
case with a trapping layer base of 800 m. The ray trace in this case is shown in Fig. 63 
with a transmitter height of 800 m. The pattern is very similar to the one observed in 
Fig. 59. One difference is that the narrow beam of focused rays slopes upward as it 
follows the path of the duct created by the trapping layer. The lower boundary of the 
ray pattern was also altered by the slope of the duct and caused the lower rays to fail to 
reach and reflect off the surface. When the transmitter was placed at 900 m and at 700 
m, the resulting ray traces showed virtually no difference in the ray patterns as compared 
to the homogeneous case. The ray patterns for a transmitter at 900 and 700 m for this 
case are shown in Figs. 64 and 65. 

3. Refractive Structure Three 

In this scenario the refractive structure was altered to simulate the presence of 

two dissimilar air masses in the region. A trapping layer at 800 m is followed 50 km 


down range by a 50 m thick trapping layer based at 900 m. The third refractivity profile 
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Fig. 59. Ray trace for homogeneous refractive structure with trapping layer and 
transmitter at 800 m. 


is a super-refractive profile with a base at 1000 m and a thickness of 100 m. The fourth 
profile is a normal refractivity proifile located 150 km down range from the original 
profile. The situation simulates that of a low pressure center slowly moving into an area 


previously dominated by high pressure where the trapping layers have formed. The low 
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Fig. 60. Ray trace for homogeneous refractive structure with trapping layer at 800 
m and transmitter at 900 in. 


pressure is simulated by the normal refractive structure observed during conditions of 
instability. 

Fig. 66 shows the effect of this varving refractive structure when the transmitter 
is positioned within the duct. The narrow bundle of rays normally associated with the 


duct appeared to begin to spread out as it passed into the super-refractive layer and then 
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Fig. 61. Ray trace for homogeneous refractive structure with trapping laver at 800 
m and transmitter at 1200 m. 


into the region of normal refractive conditions. The changing refractive structure caused 
the initial ducting to begin to Weaken with distance. This example helps to illustrate the 
focusing effect that a trapping layer has. The ray bundle is still fairly distinct 100 km 
after the trapping layer ceased to exist. A radar hole was also generated by this 


refractive profile duc to the initial trapping layer which existed out past 50km. When 
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Fig. 62. Ray trace for homogeneous refractive structure with trapping layer at S00 
m and transmitter at 700 m. 


the transmitter was placed at altitudes above and below the trapping layer the resulting 
ray patterns were once again almost identical to those shown in Figs. 60 and 62. 
4. Refractive Structure Four 
This refractive structure was identical to that described in the previous section. 


In this case the transmitter was moved so that it would be within the normal refractive 
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Fig. 63. Ray trace for upward sloping trapping layer with transmitter at 800 m. 


laver and sending the beam out towards an area where trapping layers had formed. The 
initial height of the transmitter was placed again at $00 m. The ray pattern for this sit- 
uation is illustrated in Fig. 67. The pattern is very smooth with no irregularities. The 
trapping layers do not appear to influence the pattern in any way. The 900 and 700 m 


runs produced the same type of ray trace. 
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Fig. 64. Ray trace for upward sloping trapping laver with transmitter at 900 m. 


Additional runs were made to see if the trapping layers down range would in- 
fluence any of the rays if the transmitter height was increased substantially. The pattern 
shown in Fig. 68 with a transmitter height of 1500 m shows that the ray paths are sig- 
nificantly altered by the trapping lavers resulting in a ray bundle focused downward as 


well as a large radar hole in the lower hall of thei amrent 
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Fig. 65. Ray trace for upward sloping trapping layer with transmitter at 700 m. 


5. Refractive Structure Five 
In this final case a trapping layer at $800 m with a thickness of 50 m was followed 
by a normal refractive structure 20 km away. The relatively short separation distance 
was used to represent a rapidly moving frontal system moving into an area dominated 
by high pressure. This would result in a more abrupt change from a stable situation to 


that of of an unstable deeply mixed boundary layer. Fig. 69 shows the effect of this 
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Fig. 66. Ray trace for refractive structure simulating two dissimilar air masses with 
transmitter at 800 m. 


pattern with the transmitter located in the duct at 800 m. A single narrow bundle of rays 
is formed by the trapping layer. The bundle remains relatively intact but is refracted 
along with the rest of the rays along the ray path created by the normal refractive 
Structure. One point of interest is that in this case the trapping layer did not extend far 


enough horizontally to project the rav bundle forward or to develop a radar hole. 
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Fig. 67. Ray trace for refractive structure simulating two dissimilar air masses with 
transmitter at 800 m. 


6. Significant Observations 
The refractive structure examined in this section helped to show how much of 
an effect trapping layers with a thickness of as little as 50 m can have on EM wave 
propagation. As the height of the transmitter was changed the ray paths were altered 


significantly, especially right around the trapping layer. The inhomogeneity of the 
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Fig. 68. Ray trace for refractive structure simulating two dissimilar air masses with 
transmitter at 1500 m. 


atmosphere caused some variations in the structure of the ray patterns as the rays moved 
from one refractive condition to another. The level at which trapping layers were placed 
caused an appropriate variation in the slopes of the ray bundles associated with the 
lavers. The last refractive structure helped to show that the distance over which a trap- 


ping laver exists is also essential in determining whether a radar hole will form. This is 
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Fig. 69. Ray trace for refractive structure simulating two dissimilar air masses with 
transmitter at 800 im. 


of added interest since the ranges over which EM ducting was observed during both 
MIZENX-84 and MIZEX-87 were thought to be quite short due to the many instances 


of only one ship detecting a duct. 


103 


VII. CONCLUSIONS 


This thesis examined the temporal and spatial variability of the MABL around the 
Greenland Sea MIZ and its effect on EM refractivity. The data source for the thesis was 
composed of surface observations and upper-air data measured from three ships during 
MIZEX-87. The ships were positioned relative to the ice edge in order to observe the 
changes in atmospheric variables when moving from a surface covered by ice to a surface 
covered by water. The upper-air data were composed of a set of 169 rawinsonde profiles 
collected from 19 March - 11 April 1987. 

A temporal study was done bv looking at variations in the refractive structure during 
five different synoptic regimes. During periods the MIZEX area was under the influence 
of high pressure systems, the frequency of ducting and super-refraction was highest. 
During periods dominated by low pressure, the number of refractive events decreased 
considerably. When frontal passages were recorded the boundary laver deepened to 
greater than 2000 m and the refractive structure was normal with no ducting episodes 
recorded. The large-scale forcing brought on by different mesoscale and synoptic sys- 
tems was determined to be the major element in determining the frequency of refractive 
events. This forcing also played a major role in determining the heights of the various 
refractive layers as well as their relative strengths and thicknesses. 

A diurnal study showed that incidences of ducting and super-refraction occurred 
more often during the midday and early evening launch periods. An explanation for this 
observation was not obvious and Is of sufficient interest to warrant further study. 

A spatial study was also performed. Results indicated that refractive events such 
as ducting and super-refractive lavers tended to slope upward in height relative to the ice 
edge when measured from over the ice to over the water. Greater heights over the MIZ 
were attributed to an increased occurrence of multiple vertical layers above the well- 
mixed laver for super-refractive layers. The possibility of statistical error due to insuffi- 
cient data was brought up as a possible explanation for the higher duct heights. 
Incidences of ducting were generally caused by a well formed capping inversion above 
the well mixed layer. More than half of the incidences of super-refraction were also 
caused by the development of the inversion. 

Another portion of the spatial study investigated the horizontal homogeneity of the 


refractive structure. Two or more ships observed ducts in only 35% of the cases when 
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ducts were detected. This shows that the lower atmosphere around the MIZ is relatively 
inhomogeneous when looking at that specific refractive structure. Understanding how 
the refractive structure varies with range is essential to the exploitation of the EM 
spectrum for tactical purposes. 

A ray trace model was used to show variations in the paths of EM waves caused by 
horizontal inhomogeneity in the refractive structure. [t was demonstrated that even a 
relatively thin trapping layer can greatly affect the normal ray path of a transmitting 
unit. The changes were most apparent when the transmitter was within a duct or when 
rays Were intersecting the duct in such a way as to be trapped within it. The ray paths 
also varied when the refractive structure changed from one extreme to the other. The 
most important aspect was determined to be the position of the transmitter with respect 
to the duct. Since this is a variable often controlled by the operator, knowledge of the 
location of EM ducts and other anomalous refractive layers can increase ones tactical 
advantage bv allowing the operator to place the transmitter at the best height for a given 
situation. 

The refractivity studies showed that ducting normally occurred well above the sur- 
face. The ray trace model revealed that ducting around the MIZ can cause significant 
variations in the transmission path of an airborne detection or communications system. 
An absence of surface ducting was noted in the MIZEX-87 data. This was attributed 
to higher and weaker inversion layers. Ducting should not be considered a major con- 
cern for shipboard radar systems operating around the MIZ during the spring. 

In comparing the early spring MIZEX-87 refractive structure to the summer case 
observed during MIZEX-84, it was noted that ducting occurred much less frequently 
during the spring. Duct heights were greater, duct thicknesses were less and duct 
strengths were weaker during MIZEX-87. This was attributed to a reduced level of 
subsidence in the spring as compared to the summer as well as an increased amount of 
mixing in the boundary layer caused by greater levels of heat flux. The strange dewpoint 
curves observed in many of the MIZEX-84 rawinsonde profiles were observed in only a 
small number of profiles in the MIZEX-87 data. This could be due to improvements in 
the design of the humidity sensor used in the rawinsondes. If the dewpoint phenomenon 
is not caused by instrument error it could be a seasonal phenomenon that occurs more 
frequently in the summer than in the spring. Additional research is necessary to solve 


this problem. Because the moisture in the atmosphere affects the refractive structure so 
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greatly, an understanding of this problem 1s essential to a complete understanding of the 
observed refractivity in the Arctic. 

The surge of interest in the Arctic as a theatre for military operations as well as a 
region of economic potential requires the investment of increased time and resources for 
research in this region. Development of a larger data base for the study of the Arctic 
refractive structure is a necessary part of this research as our military forces seek to es- 


tablish an increased presence in this region of the world. 
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